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1.  INTRODUCTION 


Organo-nitro  compounds  exhibit  complex  chemical  behavior,  and  resolution  of  the  mechanistic 
details  of  the  reactions  of  these  compounds  at  a  molecular  level  is  an  experimental  and  a  theoretical 
challenge.  The  chemistry  of  the  simplest  organo-nitro  compound,  nitromethane,  is  not  completely 
understood,  despite  several  experimental1'13  and  theoretical14*20  studies  done  on  the  system.  Even 
seemingly  direct  studies  of  the  unimolecular  decomposition  of  gas-phase  nitromethane  have  indicated 
complexities  in  the  mechanism  of  decomposition.2,3 

The  infrared  multiphoton  dissociation  (IRMPD)  experiments  of  Wodtke,  Hintsa,  and  Lee2  indicate 
that  isomerization  of  nitromethane  to  methyl  nitrite  competes  with  scission  of  the  C-N  bond. 

Assuming  an  isomerization  channel  and  using  RRKM  theory  to  fit  the  experimental  results,  Wodtke  et  al. 
predict  a  55.5  kcal/mol  barrier  to  isomerization  which  is  slightly  lower  than  the  C-N  bond  dissociation 
energy  (59.5  kcal/mol).  They  conclude  that  unimolecular  decay  of  nitromethane  can  occur  in  either  of 
two  competing  reactions: 

[I]  CH3N02  ->  CHJfOl  ->  CH3  +  N02  ; 

[HI  CHJJOl  ->  CH3ONO f  ->  CH30  +  NO. 

They  also  determined  the  branching  ratio  of  Reaction  II  to  Reaction  I  to  be  0.6. 

Determination  of  the  location  and  energy  of  the  transition  state  for  isomerization  of  nitromethane 
to  methyl  nitrite  has  been  a  focus  of  several  ab  initio  calculations;14*19  however,  the  multi- 
configurational  nature  of  the  ground  state  of  nitromethane  is  a  formidable  problem  in  electronic 
structure  studies.  Good  quality  ab  initio  calculations  performed  on  this  system  predict  that  the 
transition  state  structure  is  "tight"  (the  C-0  and  C-N  bond  lengths  are  approximately  2.0  A),  and  the 
energy  of  this  structure  is  approximately  76  kcal/mol  above  that  of  equilibrium  nitromethane.14,13,17  It 
has  been  suggested  that  another  transition  state  exists  that  is  more  of  a  van  der  Waals  complex,  with 
longer  C-N  and  C-0  bonds  (approximately  4.0  A).13,19  The  exact  nature  of  the  transition  state  can  be 
determined  only  by  performing  large  active  space  CASSCF  plus  Cl  calculations  using  extended  basis 
sets.  Such  multi-configurational  quantum  mechanical  studies  are  often  limited  to  systems  with  only  a 
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few  degrees  of  freedom  because  the  computational  requirements  for  these  problems  are  immense  and 
existing  computer  resources  are  limited. 

In  addition  to  the  problems  associated  with  calculating  points  on  the  potential  energy  surface, 
quantum  mechanical  dynamics  calculations  on  a  system  such  as  this  are  not  computationally  feasible. 
Classical  mechanical  methods  do  not  require  die  computer  resources  that  quantum  mechanical 
calculations  do;  therefore,  classical  trajectory  or  molecular  mechanics  calculations  are  often  employed 
to  study  the  dynamics  of  systems  with  many  degrees  of  freedom.  We  present  the  results  of  a  classical 
dynamics  calculation  of  the  unimolecular  decomposition  of  nitromethane  in  order  to  further  understand 
the  dynamics  of  a  system  with  many  degrees  of  freedom  and  complex  reaction  behavior. 

We  have  constructed  three  model  potential  energy  surfaces  to  describe  unimolecular  processes  in 
nitromethane  and  have  performed  a  series  of  classical  trajectory  studies  of  the  unimolecular 
decomposition  of  nitromethane.  Because  the  energy  and  the  geometry  of  the  transition  state  is  not 
established,  we  have  attempted  to  determine  die  effect  of  different  barrier  heights  on  the  reaction 
scheme  suggested  by  Wodtke  et  al.2  An  interesting  result  of  this  dynamics  study  is  that  all  three 
surfaces  predict  that  the  unimolecular  decay  of  nitromethane  occurs  not  just  through  Reaction  I  and  II 
but  also  by  the  reaction 

[HI]  CH^NOl  ->  CH3ONO'  CH3  +  NOr 

This  reaction  has  a  strong  energy  dependence  and  occurs  more  readily  at  higher  total  energies.  At  all 
energies,  it  makes  a  significant  contribution  to  the  number  of  isomerizing  trajectories.  The  branching 
ratio  calculated  by  Wodtke  et  al.2  was  determined  from  the  number  density  of  NO  and  N02  appearing 
as  NO+  in  the  time-of-flight  spectra.  We  assume  that  the  experiment  did  not  distinguish  between  the 
N02  formed  in  Reaction  I  from  that  which  might  be  formed  in  Reaction  m.  In  order  to  compare  our 
results  with  the  experiment,  we  define  the  branching  ratio  to  be  the  ratio  of  the  number  of  NO  formed 
to  that  of  the  number  of  N02  formed,  i.e..  Reaction  II/(Reaction  I  +  Reaction  ni). 

The  first  potential  energy  surface  (PES  1)  has  a  veiy  high  barrier  to  isomerization  (216.4  kcal/mol) 
whereas  C-N  bond  rupture  can  occur  with  60  kcal/mol.  One  might  predict  intuitively  that 
isomerization  to  methyl  nitrite  would  not  occur  on  this  surface.  However,  the  results  of  trajectories 
indicate  that  isomerization  can  occur  on  this  surface  through  a  mechanism  of  "dissociation"  to 
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CHj  +  NOj  followed  by  recombination  of  these  particles  to  form  methyl  nitrite.  Because  an  individual 
trajectory  represents  a  single  isolated  molecule  in  a  vacuum,  the  calculations  indicate  that  this 
mechanism  is  feasible  even  in  a  ” collision-free"  environment  However,  isomerization  on  this  surface 
is  rare  (between  1  to  40  isomerizations  in  every  1,000  reactive  trajectories  over  the  range  of  energies 
studied)  and  does  not  reproduce  the  experimentally  determined  branching  ratio.2  The  branching  ratios, 
as  defined  above,  are  between  0.001  and  0.025  over  the  energies  studied  for  this  potential  energy 
surface  (69.2  to  161.4  kcal/mol  above  the  zero  point  energy  of  nitromethane). 

The  second  potential  energy  surface  (PES  2)  has  an  isomerization  barrier  height  of  55.1  kcal/mol, 
the  value  suggested  by  Wodtke  et  al.2  The  third  potential  energy  surface  (PES  3)  has  an  isomerization 
barrier  height  of  47.6  kcal/mol,  the  value  predicted  from  the  MINDO/3  calculations  of  Dewar  and 
Ritchie.16  The  isomerization  of  nitromethane  to  methyl  nitrite  on  these  surfaces  does  not  occur 
through  a  dissociation-recombination  mechanism  but  rather  occurs  through  a  concerted  C-N  bond¬ 
breaking  C-0  bond-forming  process.  Although  PES  2  and  PES  3  differ  by  only  a  few  kcalAnol  in  the 
barrier  heights,  the  dynamics  results  are  dramatically  different  The  branching  ratio,  as  defined  above, 
for  PES  2  at  the  threshold  energy  of  C-N  bond  dissociation  is  0.032,  whereas  the  branching  ratio  for 
PES  3  at  the  same  energy  is  0.79.  The  branching  ratios  for  PES  2  increase  slightly  with  increasing 
energy  to  a  maximum  value  of  0.065.  whereas  the  branching  ratios  for  PES  3  decrease  with  increasing 
energy  to  a  minimum  value  of  0.12.  Additionally,  the  probability  of  Reaction  HI  increases  with 
increasing  energy  for  both  PES  2  and  PES  3,  but  the  probability  of  Reaction  n  decreases  with 
increasing  energy  for  PES  3.  Conversely,  the  probability  of  Reaction  II  is  approximately  constant  as  a 
function  of  energy  for  PES  2.  These  potential  energy  surfaces  properly  describe  the  fundamental 
vibrational  frequencies  and  geometries  for  reactants  and  products  and  the  correct  thermochemistry  for 
the  reaction  channels  shown  above.  The  main  differences  in  the  surfaces  are  in  the  transition  state 
region. 

In  addition  to  providing  mechanistic  details  of  the  unimolecular  decomposition  of  nitromethane, 
we  present  a  new  approach  to  construction  of  an  empirical  potential  energy  surface  whidh  describes  a 
many-body  system  which  can  undergo  more  than  one  type  of  reaction. 
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2.  POTENTIAL  ENERGY  SURFACE 


A  limitation  of  a  classical  trajectory  calculation,  aside  from  the  inability  to  describe  important 
quantum  mechanical  effects,  is  that  an  approximate  potential  energy  surface  must  be  used.  A  complete 
and  accurate  multidimensional  potential  energy  surface  can  be  constructed  only  by  ab  initio  quantum 
mechanical  calculations;  however,  they  are  not  computationally  possible  at  this  time  for  large 
molecular  systems.  Therefore,  the  potential  energy  surfaces  used  in  classical  mechanics  calculations 
are  usually  empirical  in  nature  and  are  based  only  in  part  on  theoretical  calculations.  Most  of  the 
potential  energy  surfaces  used  in  classical  trajectory  simulations  have  an  analytic  form  (to  allow  rapid 
calculation  of  the  forces),  with  the  parameters  of  the  potential  function  fit  to  spectroscopic,  kinetic,  and 
thermodynamic  data.  Additionally,  selected  points  calculated  by  ab  initio  methods  are  sometimes 
incorporated  into  the  fit  Usually,  most  of  the  experimental  data  correspond  to  the  equilibrium  state  of 
the  reactant  or  product  species  and  the  transition  state. 

In  our  approach  to  constructing  a  potential  energy  surface,  we  have  assumed  that  the  potential 
energy  is  the  sum  of  terms  which  correctly  describe  the  equilibrium  properties  of  each  stable  species  in 
the  system,  weighted  by  functions  which  are  dependent  on  geometric  considerations.  The  form  of  the 
overall  potential  is: 

<» 


where  i  denotes  all  the  reactant  and  decomposition  products.  The  V,  describes  the  stable  species  i  at 
equilibrium  and  is  made  up  of  terms  which  describe  stretching  and  bending  motions.  The  parameters 
of  the  functions  in  V,  are  fit  to  toe  known  data  of  the  species  near  equilibrium.  The  S,  is  a  weighting 
function  which  is  dependent  on  critical  geometries  of  the  species.  This  function  increases  or  decreases 
the  contribution  of  the  V4  to  the  overall  potential  depending  on  the  molecular  geometry. 

An  advantage  to  using  this  form  of  potential  is  that  the  weighting  functions  can  be  monitored 
throughout  a  trajectory  to  characterize  the  system  at  any  time.  As  will  be  illustrated  in  this  paper, 
following  the  weighting  functions  in  time  provides  an  easy  way  to  determine  a  reaction  mechanism  in 
a  complex  system  without  having  to  rely  on  plots  of  individual  internal  coordinates  as  functions  of 
time.  Another  advantage  to  using  this  type  of  potential  is  that  it  is  "modular;  the  different  V, 
developed  can  be  used  in  other  dynamics  studies  which  use  a  potential  energy  surface  of  the  form  in 
Equation  1.  If  other  channels  need  to  be  added,  all  that  is  necessary  is  that  the  new  V|S,  be  added  to 
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Equation  1.  and  the  existing  weighting  functions  modified  accordingly.  We  attempted  to  adjust  the 
surface  features  by  modifying  the  S,  only  and  leaving  the  V,  unchanged,  thus  developing  V,  potential 
terms  that  can  be  used  in  studies  of  other  related  systems. 

Using  the  form  of  the  potential  energy  surface  given  in  Equation  1,  the  total  potential  energy 
function  to  describe  the  nitromethane/methyl  nitrite  system  is: 

^Taul  m^CH,NOt  ^CHtNO,+^CH,ONO'  ^CH, 0W0'+^CH,0 W  $ CHtOVO 

*^CH,*NO,  ^CH,*NO,+^CH,0  *NO’  ^CH,0*N0'*^CH,0’*N0  ^CH,0‘*N<r  (2) 

Note  that  there  are  two  terms  to  describe  the  methyl  nitrite  species  and  two  terms  to  describe  the 
product  fragments  CH30  and  NO.  The  terms  differ  only  in  the  labeling  of  the  oxygens,  which  are 
labelled  O  and  O'.  These  arc  chemically  equivalent  species  and  in  reality  cannot  be  distinguished.  A 
dynamics  study  must  allow  unbiased  formation  of  methyl  nitrite  and  the  product  fragments  with  either 
O  or  O'  bound  to  the  carbon  atom.  Therefore,  it  is  necessary,  according  to  the  formulation  of  the 
potential  energy  surface  given  in  Equation  1,  to  have  terms  which  explicitly  describe  the  species  in 
which  either  O  or  O'  is  bound  to  the  carbon  atom. 

The  functions  which  make  up  the  V,  are: 

(a)  morse  functions  to  describe  bond  stretching  motion, 

V~ *  =  D,{ l-expi-aJR-RJ)}2  -  D,;  (3) 

(b)  harmonic  oscillator  functions  to  describe  bending  motion, 

V**  =  #9  -  8.  /  ;  (4) 

(c)  and  truncated  cosine  series  to  describe  torsional  motion,  where  applicable, 

Vj(x)  =  I  A,cos(nx).  (5) 

Where  necessary,  these  standard  functions  have  been  modified  as  described  below. 
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Because  the  S,  is  a  weighting  function  for  the  contribution  of  the  V,  to  the  overall  potential,  it 
must  be  a  smoothly  varying  function  between  0  and  1.  Additionally,  its  first  derivatives  must  be 
continuous.  Although  there  is  substantial  information  about  die  geometry  and  intramolecular 
interactions  of  each  reactant  or  product  species  near  its  equilibrium,  little,  if  anything,  is  known  about 
the  regions  of  the  potential  between  the  limits  of  the  reactant  and  product  species.  Therefore,  we  have 
little  information  to  guide  us  in  choosing  the  functional  form  of  the  S,  to  describe  the  regions  between 
the  asymptotic  limits  of  reactants  and  products. 

The  function  which  makes  up  the  S,  is: 

F'(Q)  =  exp[-G(Q-Qj),  (6) 

where  Q  defines  some  internal  coordinate  and  N  is  an  even  integer  greater  than  or  equal  to  2. 

Figures  la  and  lb  show  F'  for  different  N  and  G,  respectively.  A  modified  form  of  F'  used  in 
constructing  the  S,  is: 


F(Q )- 


Q2Qn 


(7) 


and  is  illustrated  in  Figure  lc. 

The  description  of  the  details  of  VXaMl  will  be  as  follows:  each  reactant  or  product  fragment  will 
be  discussed  in  separate  sections,  and  the  V,  and  S,  will  be  detailed,  along  with  the  experimental  data 
and  theoretical  information  about  the  fragment  which  are  used  in  fitting  the  parameters.  The  functions 
F'  and  F  will  be  denoted  with  subscripts,  for  use  in  referencing  the  parameters  given  in  Table  1. 

2.1  Nitromethane. 

2.1.1  Form  of  the  Potential.  The  term  which  describes  equilibrium  nitromethane  is: 

V  a  E  v *'***+  y Jfcrte*  y  ytaMc* 

VCH,NO,  **  VCH  *  VCN  *  **  V NO 

*  £  Vhch  *  Vowo  +  [  E  +  E  Van  ]  F/R^,  (8) 
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where  all  of  the  functions  have  been  defined  in  Equations  3-7.  The  form  of  the  function  which 
deteimines  the  contribution  of  VCH3N02  to  the  overall  potential  is  a  function  of  Rq,  and  the  two  C-N-0 
angles,  ®CNO  and  0cno-: 

& CH,NO,m  CNO^ «(®cwo  )  •  @) 

This  fiinction  attenuates  VCH3NC>2  as  the  molecule  approaches  decomposition  to  CH3  +  NOj  or 
isomerization  to  either  CH3ONO '  or  CHjO'NO.  The  parameters  for  the  Fj  functions  are  given  in  Table  1. 

2.1.2  Experimental  Data  and  Theoretical  Information  Used  for  Parameter  Fitting. 

(1)  Structure.  Gas  phase  electron  diffraction  and  microwave  spectra  of  CHjN02  have  been  used 
by  several  groups44  to  determine  the  structural  parameters  of  CH3N02  and  to  study  the  internal  rotations 
of  the  methyl  group  relative  to  the  N02  group  in  nitromethane.  While  the  interest  of  die  early  microwave 
studies  focused  on  determination  of  the  barrier  height  to  internal  rotation,44  Cox  and  Waring7  provide  an 
experimental  structure  for  nitromethane  (see  Table  2).  The  analysis  of  the  experimental  results  is  based 
on  assumptions  about  the  methyl  top.  Although  the  heavy  atom  bond  lengths  were  accurately  determined 
in  this  experiment,  the  C-H  bond  length  was  not  determined  but  was  assumed  to  be  1.089  A.  A  later 
microwave  spectral  study  of  CHjN02  was  performed  by  Cox  and  co-workers*  to  determine  the  structure 
of  the  methyl  group.  The  study  shows  that  Rc„  =  1.0979  ±  0.0009  A,  and  0^  =  107.61  ±  0.15°. 

Electronic  structure  calculations14,13  show  a  dependence  of  geometric  parameters  on  orientation  of 
the  methyl  group  with  respect  to  the  N02  group.  Changes  in  bond  lengths  and  the  degree  to  which  the 
methyl  and  N02  rotors  tilt  off  die  C-N  bond  axis  as  the  nitromethane  molecule  undergoes  internal  rotation 
are  predicted  by  the  ab  initio  calculations.  McKean  and  Watt,9  in  an  interpretation  of  infrared  (IR) 
spectra  of  nitromethane,  predict  changes  in  the  C-H  bond  lengths  of  up  to  0.006  A  during  internal 
rotation,  affirming  the  effects  predicted  by  first  principles  calculations.  The  barrier  to  rotation  of  the 
methyl  group  about  the  C-N  bond  axis  is  0.006  kcal/mol,6  causing  it  to  be  a  free  rotor  at  room 
temperature.  Because  we  are  specifically  interested  in  high  energy  processes  in  nitromethane 
(unimolecular  decomposition)  and  the  geometric  changes  due  to  the  rotation  are  small,  we  have  chosen 
not  to  include  terms  in  Equation  8  which  are  dependent  on  internal  rotation.  These  effects,  however, 
could  easily  be  incorporated  into  the  potential  by  using  the  form  of  the  potential  suggested  by  McKean 
and  Watt,9  which  is  a  function  of  torsional  angle. 
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Table  2.  Structural  Parameters  of  Nitromethane 


Internal 

Coordinate* 

ExptT.b 

PES  le 

PES  2C 

PES  3C 

R(CH) 

1.088 

1.076 

1.076 

1.076 

R(CN) 

1.489 

1.479 

1.479 

1.479 

R(NO) 

1.224 

1.221 

1.221 

1.221 

Z(HCN) 

107.2 

107.40 

107.40 

107.40 

Z(ONO) 

125.3 

125.80 

125.87 

125.87 

*  Units  are:  R  in  A,  Z  in  degrees 
b  Cox  and  Waring,  Ref.  7 

*  G.  F.  Adams,  unpublished  results. 


(2)  Vibrational  Analyses.  Nitromethane  has  14  fundamental  vibrational  modes  and  one  internal 
rotational  degree  of  freedom  (torsional  mode).  Assignments  of  frequencies  from  infrared  (IR)  spectra  are 
provided  by  several  groups9-13  and  shown  in  Table  3. 

(3)  Bond  Dissociation  Energies.  The  dissociation  energy  t-f  the  C-N  bond  in  nitromethane  has 
been  measured  to  be  58  to  60  kcal/mol,21  while  the  C-H  and  N-0  bond  dissociation  energies  are  not 
firmly  established.  The  C-H  bond  energy  in  nitromethane  has  been  measured  at  107.5  kcal/mol,22  but 
McKee13  predicts  a  substantially  lower  bond  dissociation  energy  of  88  kcal/mol  based  on 
MP2/6-31G*  +  ZPC  level  and  MCSCF/6-31G  +  Cl  level  calculations.  We  have  assumed  a  C-H  bond 
dissociation  energy  of  102  kcal/mol  as  a  reasonable  approximation.  The  N-0  bond  dissociation  energy  in 
CHjN02  has  been  estimated  to  be  93.5  kcal/mol.23 

2.1.3  Results  of  the  Fit  The  parameters  of  the  potential  energy  function  described  in 
Equations  8-9  were  fitted  to  the  geometric  parameters  predicted  by  SCF  calculations  of  Adams,14  except 
for  the  slight  differences  due  to  internal  rotation.  His  calculated  results  are  in  good  agreement  with  the 
experiment,  except  for  the  C-N  bond  distances.  The  internal  coordinates  of  CH3N02  at  equilibrium 
predicted  by  the  three  potential  energy  surfaces  are  listed  in  Table  2. 


9 


Table  3.  Normal  Mode  Frequencies1  for  Nitromethane 
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Frequencies  are  in  cm'1 


Assuming  the  bond  dissociation  energies  listed  above,  all  other  parameters  in  Equations  8-9  were 
fitted  to  the  experimental  frequencies  shown  in  Table  3.  The  three  surfaces  predict  13  vibrational 
modes,  shown  for  comparison  with  experimental  values  in  Table  3.  No  torsional  mode  exists  as  this 
motion  was  not  explicitly  described  in  Equations  8-9,  as  explained  above.  Additionally,  the  N02  out- 
of-plane  rock  is  not  present  The  calculated  frequencies  are  in  reasonable  agreement  with  the 
experiment  except  for  vibrational  modes  ol0  and  o4.  The  best  fit  parameters  are  given  in  Table  4. 

2.2  Methyl  Nitrite. 


2.2.1  Form  of  the  Potential.  The  term  which  describes  a  methyl  nitrite  molecule  at  equilibrium 
is  considerably  more  complex  than  that  of  the  nitromethane  molecule.  Two  geometric  isomers  of 
methyl  nitrite  exist— The  "cis"  and  "trans"  species.  The  force  constants  and  geometries  differ  in  these 
isomers.  Also,  there  is  evidence  that  forces  and  geometries  change  upon  internal  rotation  of  the 
methyl  group  in  methyl  nitrite.  Using  the  available  information,  we  have  incorporated  the  effects  of 
internal  rotation  in  the  potential  energy  function.  The  symbol  i'  will  be  used  to  denote  the  CONO 
dihedral  angle,  and  x„  i  «  1,  2,  or  3,  will  denote  the  H0)CON  dihedral  angles.  Without  loss  of 
generality,  we  will  use  the  VCH30NOSCH30N0.  term  to  describe  the  methyl  nitrite  potential;  the 
vch30'N(>Sch30'N0  term  has  the  same  form,  except  for  tire  position  of  the  oxygen  atoms.  Additionally,  a 
few  terms  differ  between  PES  1  and  PES  2  and  3.  The  difference  occurs  only  in  use  of  attenuation 
functions  and  die  parameters  of  the  attenuation  functions.  The  differences  will  be  specified  below. 


The  forms  of  the  potential  describing  methyl  nitrite  at  equilibrium  are: 
PES  1 


V  m 

vCHtONO ' 

+ 


wsr  (t  *  >+vsr\%  o  ♦  v%r  c*  o 


,  Struck  . 


V£T  (X  0  +  ( +V„ch^0 

"d)CKm  ”<2)  c”(3> 


V—a,™  * 


*  t  ,X1,Tj,Tj)+Vcav(X  j/tjiTj) 

♦  vr  (tO]F7(*w)F5(*co). 


(10a) 
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Table  4.  Potential  Energy  Parameters,  Nitromethane  Term 


Morse  Parameters 

Bond 

D.(eV) 

a  (A1) 

R*(A) 

CH 

4.423174 

1.89 

1.076 

CN 

2.601867 

2.085 

1.479 

NO 

4.045904 

2.377 

1.221 

Harmonic  Oscillator  Parameters 

Angle 

k  (eV/rad2) 

0«,  (degrees) 

HCH 

1.375 

111.46 

HCN 

2.45 

107.40 

CNO 

3.20 

117.10 

ONO 

9.35 

125.80 

-Yvsra  '.x  o+v^cx  o 

.  ,  T  Struck  x/Bmd  >  .  *  jBtnd  /,* 


+VjST( X')*KMch1x'^  +  K  ch^ 

"(1)  CW(2)  "(2)  CW(3) 


+VH  CHX'^  +  <  E  vr*jtw 

"(I)**®  w(0cd 

+  t^iawo*  (x  ,x,t  x2,Tj)  +  Vcawx  .XpXj,  x3) 

♦  VI<xOJF7(i?aw)}P’s(J?co). 
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We  have  indicated  the  dependence  of  each  term  on  the  torsion  angles.  These  terms  are,  of  course, 
dependent  on  other  internal  coordinates.  Those  dependencies  are  not  explicitly  shown.  Note  that  the  only 
difference  between  Equation  10a  and  Equation  10b  is  that  the  HCH  bending  terms  described  in  Equation  10a 
are  not  attenuated  by  the  function  FjO^o). 

The  modified  Morse  function  which  describes  the  C-H  stretching  motion  is: 

vStftr.x)  -Da,(1.0-exp(-aal[*OT(i)-*“l)F-Do,.  (U) 


where  da,  is  a  function  of  t',and  s;  and  has  the  form: 

actl  -  A1ch+B1c"+A2cw[1.0-cos(t')] 

+  {  Cf+C™ [1.0-cos(t  0  ] } cos2^ ,). 

The  modified  Morse  function  which  describes  the  C-0  stretching  motion  is 

-  Dco  {1.0  -exp(-a CO[RCO-R™]Y-DCO, 

where  is  a  function  of  x', 

aco  -  B?°+B2co[  1.0-cos(x')]. 

The  O-N  bond  is  described  by  a  modified  Morse  function, 

v£T*(t  0  -  ft.0  -  exp  (-<xw  [Rw-R?} )  P 

where  is  a  function  of  x', 

aoN  -  *1w+02oiv[1.O-cos(t')]. 


(12) 


(13) 


(14) 


(15) 


(16) 
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The  terminal  N=0  bond  is  described  by  a  modified  Morse  function, 

v£T(t  0  -  DNw0  (1.0-exp(-c *  (17) 

where  is  a  fimction  of  x', 

aNwt}  -  Bf^+fl2^[lX)-a>sCt')].  <18> 

The  potential  term  describing  motion  of  H(1)CHP)  angle  is 

cHpj-®  J  *  (19) 

where  0O  is  a  function  of  internal  rotations  and  has  the  form 

0O  -  rf^^^l.O-cosfx  0] 

+  {B1hch+B?ch  [  1  .0-cos(x  01  fcos  (Xj) 

+  (C/^+C/^ll.O-cosfx  01  )cos2(x2).  (20) 

The  terms  describing  the  motion  of  the  H^CH^  and  the  H^CHy,  angles,  respectively,  have  the  same  form 
as  Equations  19  and  20,  except  for  0O  for  the  H^CH^  angle  is  a  fimction  of  x,  and  0O  for  the  H(1)CH0) 
angle  is  a  fimction  of  tj. 

The  ONO  bending  potential  term  has  the  form 

Vomft  tX{lvX^  *  ^OWo[®OWO~® J2’  (71) 

where 

*«vo  -  *r°+  C^t  1.0-cos  (x  }] .  (72) 
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and 


ef®  -  rr°^ni.o-cos(to] 

3 

+te,aw[1.0-cos(T  o  ] }  n  [1.0-cos(T  )].  (23) 

i-1 

The  CON  bending  potential  has  the  form 

•'Cj.'tj.Xj)  -  (24) 

where 

*cav  «  ^C<W+A,CW[1.0-COS(t  0]  ,  (25) 

and 

3  3 

e,  -  7’/**'  +  (  b™  n  [Lo-cosctpi+o^n  [l.o+cosCx,.)]} 
i-1  i-1 

3 

X  [  1.0-cosCO]  +  lCyCON  n  (l.O+cosd^lHl.O  +  cos(x0  ].  (26) 

i-1  . 

The  potential  term  which  describes  the  H^CO  angle  bending  potential  is 

VH^co  (X  '.X .)  -  *W(,)CO  [0W(j)co  'e  J*  •  (27) 

where 

*«(0c<>  •  1.0-cos  (x  oi 

+  ( cr+C™  [1.0-cos(x  o  ] }  cos(x  )  ,  (28) 
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and 


9,  -  7'*WC0+  A^°*  A^°  [1.0-cos(x  OJ 
+  to?00  +  B™  [1.0-cos(t')]}cos(xj) 

+  {C™  *  Cf*  [1 .0-cos(x  0 1 }  cos2  (x  ()  .  (29) 

The  potential  which  describes  the  internal  rotation  of  the  CONO  group  is  that  given  by  Darsey  and 
Thompson24  and  has  the  form  given  in  Equation  5. 

The  weighting  function  SCHj0S0-  is  a  function  of  Rcc  Ron1  Rco  -  and  0^  and  has  die  following 
forms  for  the  three  surfaces: 

PES  1 

Sat 3  aw-  «  ^co)  W^'9(U  [l.O-F^J?*,.)].  (30a) 

PES  2  and  3 

^CHj ono'  -  F2(/?c^F3(Fcw)[1.0-F|(eaw)][1.0-F1#(/?co.)]  .  (30b) 

This  function  attenuates  VCH30N0.  as  the  molecule  approaches  decomposition  to  CH,  +  N02  or  CH30  + 
NO',  or  when  the  system  is  in  a  geometry  consistent  with  nitromethane.  The  weighting  function  S^o^o 
has  a  similar  form  as  Equation  30: 

PES  1 

^CHj  OVO  "  ^2 (*CO-) ^3 (*0 tl.0-F16  (Fco)]  .  (31a) 

PES  2  AND  3 

ScHj  o*o  "  F2(Fco.)F,(Fow)[l.O-Fl(0cw.)]  [1.0-F16(Fco)1  .  (31b) 
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Parameters  of  the  weighting  functions  Ft  are  given  in  Table  1. 

2.2.2  Experimental  Data  and  Theoretical  Information  Used  for  Parameter  Fitting. 

(1)  Structure.  Accurate  structures  of  both  cis  and  trans  isomers  of  methyl  nitrite  have  been 
derived  from  microwave  spectral  studies23,2*  as  well  as  barrier  heights  to  internal  rotation.23,26 

The  experimental  results  show  decreases  in  the  CON  and  ONO  bond  angles  upon  cis-trans  isomerization 
as  well  as  changes  in  the  bond  lengths  and  the  tilt  of  the  methyl  group.  The  methyl  tilt  is  defined  as  the 
angle  between  the  C-0  bond  and  the  perpendicular  from  die  carbon  atom  to  die  basal  plane  of  the  CH3 
pyramid.  Turner  et  al.23  determined  the  methyl  tilt  to  be  5.4°  ±1.5°  away  from  the  nitrosyl  oxygen  atom 
in  the  cis  isomer,  the  structure  and  conformation  of  the  methyl  group  for  trans-CH3ONO  was  not 
determined.  The  decreases  in  the  CON  and  ONO  bond  angles  upon  conversion  from  the  cis  to  trans 
isomer  is  thought  to  be  due  to  reduced  steric  interaction  in  the  trans  conformation.  The  experimentally 
determined  structural  parameters  of  both  cis-  and  trans-methyl  nitrite  are  given  in  Table  5.  Ab  initio 
calculations  of  methyl  nitrite  predict  changes  in  bond  lengths  and  valence  angles  upon  internal 
rotations.13,27  McKee13  provides  calculated  changes  in  the  methyl  tilt  upon  internal  rotation  and  discusses 
these  in  the  context  of  conformational  stabilities  within  methyl  nitrite.  Ha,  Meyer,  Ghosh,  Bauder,  and 
Gflnthard27  examined  the  structural  relaxation  of  the  cis  conformer  upon  the  internal  rotation  of  the  methyl 
group  and  its  effect  on  the  rotational  constants  upon  torsional  excitation.  These  changes  in  geometric 
parameters  upon  internal  rotation  are  shown  in  Table  5. 

(2)  Vibrational  Analyses.  Methyl  nitrite  has  15  vibrational  degrees  of  freedom.  The  IR  gas  phase 
spectrum  of  cis-  and  trans-methyl  nitrite  recorded  over  the  range  4,000-300  cm'1  is  reported  by  Ghosh  and 
GOnthard,2*  and  they  provide  assignments  of  the  vibrational  frequencies  for  both  the  cis  and  trans  isomers. 
These  assignments  are  given  in  Tables  6  and  7. 

(3)  Bond  Dissociation  Energies.  Because  the  C-H  bond  dissociation  energy  is  not  well 
established,  we  have  assumed  a  value  of  102  kcal/mol,  die  same  value  assumed  for  nitromethane.  The 
O-N  bond  energy  has  been  measured  to  be  39  to  41  kcal/mol;21  we  used  41.2  kcal/mol.  Based  on  heats 
of  formation  calculated  by  Melius,17  we  assume  the  bond  dissociation  energy  of  the  terminal  N=0  bond 
is  151.8  kcal/mol  and  that  of  the  C-0  bond  is  58  kcalAnol. 
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Table  5.  Structural  Parameters  of  Cis-  and  Trans-Methyl  Nitrite 


Cis-Methyl  Nitrite  | 

Int.® 

Exptb 

Theor.® 

PES 

PES 

PES 

Theor.® 

PES 

PES 

PES 

Coord. 

1 

2 

3 

1 

2 

3 

Staggered 

Eclipsed 

0=N 

1.182 

1.162 

1.162 

1.168 

1.168 

1.160 

1.162 

1.166 

1.166 

N-0 

1.398 

1.326 

1.326 

1.348 

1.349 

1.331 

1.326 

1.338 

1.338 

u 

1 

o 

1.437 

1.421 

1.421 

1.435 

1.435 

1.425 

1.421 

1.429 

1.429 

C-H, 

1.09 

1.078 

1.080 

1.075 

fllll 

C-H0 

1.102 

1.080 

n 

S1S1 

ZONO 

114.8 

115.2 

115.2 

118.0 

118.0 

116.4 

116.4 

118.0 

118.0 

ZNOC 

IBS 

117.6 

117.6 

122.6 

122.6 

120.7 

120.7 

123.7 

123.7 

ZOCH, 

101.8 

104.4 

104.4 

104.7 

104.7 

111.7 

111.7 

111.3 

111.3 

ZOCH0 

109.9 

110.8 

110.8 

110.6 

110.6 

107.2 

107.2 

107.4 

107.4 

ZHoCH0 

108.1 

— 

108.1 

107.7 

107.7 

— 

111.6 

112.3 

112.3 

HjCH0 

— 

— 

111.3 

111.6 

111.6 

— 

109.6 

109.1 

109.1 

Trans-Methyl  Nitrite  jj 

0=N 

1.164 

1.156 

1.162 

1.160 

1.160 

1.155 

1.162 

1.160 

1.160 

N-0 

1.415 

1.339 

1.326 

1.326 

1.326 

1.339 

1.326 

1.326 

1.326 

O-C 

1.436 

1.415 

1.421 

1.421 

1.421 

1.416 

1.421 

1.421 

1.421 

C-H, 

1.09 

1.078 

1.081 

El 

1.080 

C-H. 

1.09 

1.082 

1.080 

1.080 

Ed 

1.080 

ZONO 

111.8 

111.5 

111.5 

111.5 

111.5 

111.6 

111.5 

111.5 

111.5 

ZNOC 

109.9 

111.3 

111.3 

111.3 

111.3 

112.6 

112.6 

112.6 

112.6 

ZOCH, 

109.5 

105.6 

105.6 

105.6 

105.6 

110.1 

110.1 

110.1 

110.1 

ZOCH. 

109.5 

110.7 

110.7 

110.7 

110.7 

108.4 

108.4 

108.4 

108.4 

ZH.CH. 

109.5 

— 

108.2 

108.2 

108.2 

— 

110.5 

110.5 

110.5 

ZH.CH. 

109.5 

— 

110.8 

110.8 

110.8 

— 

109.7 

109.7 

109.7 

*  Unit*  are:  R  in  A,  Z  in  degree*.  H,  designate*  the  hydrogen  atom  in  the  plane  of  the  molecule  and  H.  designates  the 
hydrogen  atom  out  of  the  plane  of  the  molecule. 

b  Turner,  CorkilL,  and  Cox,  Reference  25. 

*  McKee,  Reference  15. 
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Table  6.  Normal  Mode  Frequencies1  for  Cis-Methyl  Nitrite 


Mode 

Expt.b 

PES  1 

PES  2 

PES  3 

PES  1 

PES  2 

PES  3  | 

Staggered 

Eclipsed 

3038 

3051.7 

3069.1 

3069.1 

3044.9 

3051.1 

3051.1 

2993 

3020.2 

3018.8 

3018.8 

3027.2 

3027.9 

3027.9 

« i 

2956 

2913.9 

2914.7 

2914.7 

2912.8 

2914.7 

2914.7 

1620 

1616.4 

1683.9 

1682.6 

1616.5 

1621.1 

1620.8 

1458.5 

1472.8 

1550.2 

1549.0 

1483.4 

1534.5 

1533.9 

1444.7 

1437.5 

1444.9 

1444.4 

1453.7 

1451.5 

1451.5 

1410 

1408.1 

1393.2 

1393.2 

1441.5 

1448.1 

1448.1 

1240 

1270.3 

1236.2 

1236.1 

1253.2 

1243.4 

1242.9 

«13 

1003.5 

9752 

1182.1 

1179.7 

1027.3 

1040.4 

1040.0 

*>7 

990.5 

937.9 

956.4 

955.8 

849.2 

977.6 

977.5 

841 

708.0 

878.1 

877.0 

756.5 

843.6 

843.3 

628 

629.9 

653.5 

653.3 

631.0 

687.3 

687.3 

V10 

346.2 

356.3 

619.1 

619.0 

375.3 

636.6 

636.6 

t  (CONO) 

— 

347.3 

313.3 

313.1 

353.0 

376.1 

376.0 

T  (HCON) 

— 

— 

176.0  i 

175.9  i 

— 

126.7 

126.7 

*  Frequencies  an  in  cm1. 
b  Ghosh  and  GOnthaid,  Ref.  28 
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Table  7.  Nonnal  Mode  Frequencies1  for  Trans-Methyl  Nitrite 


*  Frequencies  are  in  cm1. 

*  Ghoeh  and  GOnthard,  Ref.  28. 
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2.2.3  Results  of  the  Fit  Because  McKee13  provided  the  most  detailed  information  about 
geometric  changes  in  methyl  nitrite  upon  internal  rotation,  we  have  fitted  the  parameters  of  the  functions 
in  Equations  10-31  to  his  ab  initio  values. 

Using  the  bond  dissociation  eneigies  discussed  above,  all  other  parameters  were  fitted  to  give 
agreement  with  the  experimental  vibrational  frequencies  shown  in  Tables  6  and  7.  Vibrational  frequencies 
calculated  from  a  normal  mode  analysis  using  these  three  surfaces  are  provided  for  comparison  with 
experiment  in  Tables  6  and  7. 

PES  1  predicts  14  vibrational  modes,  including  the  torsional  mode  of  the  CONO  dihedral  angle. 
The  frequency  of  this  motion  has  not  been  observed  experimentally  for  either  the  cis  or  trans  isomer. 
Agreement  with  experiment  is  good  with  the  exception  of  mode  u,  for  the  cis  isomer  and  mode  u13  for 
the  trans  isomer. 

PES  2  and  PES  3  predict  13  vibrational  modes  for  cis-methyl  nitrite  and  14  vibrational  modes  for 
trans-methyl  nitrite.  The  normal  frequencies  for  cis-methyl  nitrite  differ  from  those  calculated  using  PES 
It  particularly  for  modes  v4,  uJ3,  v,,  and  \j,0  for  the  cis-staggered  species,  and  o4,  u,,  u,.  and  o,0  for  the 
cis-edipsed  species.  Because  the  cis-staggered  methyl  nitrite  has  an  imaginary  frequency  for  the  torsional 
motion  of  the  methyl  group  about  the  C-0  bond,  this  indicates  it  is  a  transition  state  species  on  PES  2  and 
PES  3.  The  eclipsed  cis-methyl  nitrite,  however,  has  no  imaginary  frequencies  and  therefore  is  a  stable 
configuration  on  all  three  surfaces.  Vibrational  frequencies  of  equilibrium  trans-methyl  nitrite  are  the  same 
for  all  three  potential  energy  surfaces.  Agreement  with  experiment  is  good,  except  for  v13.  The  best  fit 
parameters  for  the  V,  are  given  in  Table  8. 

2.3  Product  Fragments  Methvl  and  Nitrogen  Dioxide  Radicals. 

2.3.1  Form  of  the  Potential. 

The  +  noj  term  is  written  as: 


/  at  V 

fCMi*NOl  yCH3 


rw2* 


(32) 
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Table  8.  Potential  Energy  Parameters,  Methyl  Nitrite  Team 


a  (A1) 

k  (eV/rad2) 

6a,  (degrees) 

A  CH 
"1 

-0.023930 

^ HCH 

1.585558 

tHCH 
*  o 

108.11 

A  CH 

a2 

-0.041760 

.OHO 

8.533768 

A  HCH 

A\ 

2.522218 

B™ 

1.870784 

cr 

0.805765 

aT 

-0.038579 

C™ 

0.023930 

u CON 
* 0 

5.355633 

bhch 

1.777910 

nCH 

'■'2 

0.007841 

A™ 

-1.081012 

b2kh 

ESB 

B?° 

2.704537 

dhco 

1.789699 

cHCH 

-0.758857 

B2 

0.033294 

a  HCO 
"1 

-1.405667 

HCH 

-0.212692 

B™ 

3.355830 

A  HCO 

a2 

0.355721 

tONO 
*  0 

116.4 

B™ 

-0.428253 

cr 

1.405667 

A  OHO 

-2.45 

0N-0 

Bl 

2.225035 

fiHCO 

C2 

-0.384331 

Br 

-1.20 

117.6 


-6.30 


-5.0 


3.10 


104.4 


4.915147 


0.401516 


3.636363 


-0.686364 


-1.260600 


-0.506062 


The  function  used  to  describe  the  methyl  radical  is: 


-  E  VST*  E  {  V™  [\.0-Fs(Rco)][\.0-Fs(Rco.)  ] },  (33) 

where  all  the  terms  are  defined  in  Equations  3-7.  Several  others  have  fit  analytic  functions  that 
reproduce  experimental  frequencies  for  CH3.29  Some  of  these  functions  have  included  terms  which 
explicitly  describe  the  out-of-plane  large  amplitude  (OPLA)  motion  of  the  methyl  radical.  Because  we 
think  addition  of  this  term  probably  will  not  affect  the  dynamics  of  the  reactions  of  CHjNOj,  we  have 
opted  to  leave  it  out  However,  a  few  groups  ™  provide  terms  which  describe  this  motion  that  could 
easily  be  incorporated  into  this  potential 

The  form  of  the  term  used  to  describe  N02  is: 


£  +  v££  [1.0-Fj  (Rce)  ]  [1.0-Fj  (Fco.) ). 


(34) 


where  all  die  terms  are  defined  in  Equations  3-7.  lire  weighting  function  for  VCHj  +  NC>2  has  the  form: 

^chj  *no2  m  t ]  [1.0  F2  (/?co)  ]  [  1.0-Fj  (Rco)  ]Fi(RN(j)Fi(RNO).  (36) 

Parameters  for  the  weighting  functions  F,  are  given  in  Table  1. 

2.3.2  Experimental  Data  and  Theoretical  Information  Used  for  Parameter  Fitting. 

(1)  Structure.  The  structure  of  the  methyl  radical  in  the  ground  state  was  established  by 
Herzberg;30  it  is  planar  with  symmetry.  Experimentally  determined  values  of  the  structural 
parameters  are  given  in  Table  9.  The  structural  parameters  of  ground  state  N02  were  determined  from 
microwave  spectra31  and  are  also  given  in  Table  9. 
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Table  9.  Structural  Parameters  of  CH,  +  N02 


Internal  Coordinate* 

Expt’l.^ 

PES  1 

PES  2 

PES  3 

Methyl  radical  j 

R(CH) 

1.079 

^7' 

ZHCH 

120.0 

120.0 

120.0 

120.0 

|  Nitrogen  Dioxide  radical  | 

R(NO) 

1.1934 

1.1934 

1.1934 

1.1934 

ZONO 

134.07 

134.4 

134.5 

134.5 

*  Units  are:  R  in  A,  Z  in  Degrees 

b  Methyl  radical:  Herzberg,  Ref.  30 

*  Nitrogen  Dioxide  radical:  Bird  et  al„  Ref.  31. 


(2)  Vibrational  Analyses.  Although  most  of  the  experimental  studies  of  the  vibrational  spectrum 
of  CH3  have  focused  on  the  OPLA  motion,32  Jacox33  has  provided  a  compendium  of  vibrational  data 
for  the  ground  state  methyl  radical,  which  are  listed  in  Table  10.  The  fundamental  vibrational 
frequencies  of  N02  have  been  determined  by  Bird  et  al.su  and  by  Sams  and  Lafferty;34  they  are  listed 
in  Table  10. 

(3)  Bond  Dissociation  Energies.  Although  the  bond  dissociation  energies  of  the  product 
fragments  have  not  been  firmly  established,  the  thermochemistry  of  the  unimolecular  dissociation  of 
nitromethane  to  CH3  +  N02  has  been  determined.33  Experiments  indicate  that  nitromethane  is 
approximately  60  kcal/mol  more  stable  than  the  product  fragments;  using  this  information  and  Melius’ 
calculated  values17  for  N-0  bond  energies  in  N02,  we  must  use  a  C-H  bond  dissociation  energy  of 
116.2  kcal/mol  to  obtain  the  correct  thermochemistry.  This  is  not  an  unreasonable  value  for  a  C-H 
bond. 

2.3.3  Results  of  the  Fit  The  internal  coordinates  at  equilibrium  of  the  methyl  and  nitrogen 
dioxide  radicals  predicted  from  the  three  potential  energy  surfaces  are  shown  in  Table  9.  Assuming 
the  bond  dissociation  energies  given  above,  all  other  parameters  were  fitted  to  give  agreement  with  the 
experimentally  observed  spectroscopic  data.  Normal  mode  analysis  of  the  product  fragments  gives  the 
fundamental  frequencies  shown  in  Table  10.  The  three  surfaces  predict  five  fundamental  frequencies 
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Table  10.  Normal  Mode  Frequencies*  for  CH3  +  N02 


|  Mode 

ExpLb 

PES  1 

PES  2 

PES  3 

]  Methyl  radical  | 

o(CH) 

3160.821 

3176.0 

3176.0 

3176.0 

o(CH) 

3004.8 

2988.6 

2988.6 

2988.6 

5(CHj) 

1396® 

1404.2 

1404.2 

1404.2 

OPLA 

606.4d 

— 

— 

— 

|  Nitrogen  Dioxide  radical  | 

Mode 

Expt* 

Exptf 

PES  1 

PES  2 

PES  3 

o(NO) 

1665.5 

1616.8 

1620.3 

1620.3 

1620.3 

o(NO) 

1357.8 

1319 

1220.3 

1220.3 

1220.2 

5(ONO) 

756.8 

749.6 

766.4 

766.4 

766.4 

Frequencies  are  in  cm.'1 

Jacox,  Ref.  33 

CH,  in  rare  gas  matrix 

Out-of-Plane  Large  Amplitude  Motion;  see  text 

Bird  et  aL,  Ref.  31 

Sams  and  Lafferty,  Ref.  34. 


for  CH3  rather  than  six  because,  as  discussed  above,  a  term  which  explicitly  describes  the  OPLA 
motion  was  left  out  Agreement  with  experiment  for  the  other  frequencies  is  good  to  within  20  cm1. 


Fundamental  frequencies  of  nitrogen  dioxide  calculated  from  these  potentials  are  also  shown  in 
Table  10.  Agreement  with  experiment  is  good  except  for  the  symmetric  stretch,  which  differs 
substantially  from  the  experimental  value.  The  calculated  equilibrium  geometry  and  fundamental 
frequencies  do  not  differ  among  the  three  potential  energy  surfaces.  The  best  fit  parameters  for  the  V, 
are  given  in  Table  11. 


2.4  Product  Fragments  Methoxv  and  Nitrogen  Oxide  Radicals. 

2.4.1  Form  of  the  Potential.  As  in  the  description  of  the  methyl  nitrite  potential,  we  choose, 
without  loss  of  generality,  to  describe  the  methoxy  and  nitrogen  oxide  potentials  using  the  term 
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Table  11.  Potential  Energy  Parameters,  CH3  +  N02  Term 


|  Morse  Parameters 

Bond 

D.(eV) 

a(A-‘) 

R*,(A) 

CH 

5.038950 

1.8125 

1.080 

NO 

3.122241 

2.90 

1.1934 

|  Harmonic  Oscillator  Parameters 

Angle 

— 

k  (eV/rad2) 

0*,  (degrees) 

HCH 

— 

1.282 

120.0 

ONO 

— 

6.25 

134.4 

VcHjO  +  NO* 


V  mV  +V 

yCHiOMO'  V  CH^O  NO' 


The  potential  describing  the  methoxy  radical  is 

tStntch  |  r  Strttck 


v  -V  tyBmd  yBtdx 

VCH^O  Z'VCH  *V CO  +**'vHCH  +  vHCO  •' 


(36) 


(37) 


where  the  terms  in  this  equation  are  defined  in  Equations  3-7.  The  term  describing  the  nitrogen  oxide 
species  is  simply: 


’NO 


y  Struck 

/NO  » 


(38) 


which  is  defined  in  Equation  3. 

2.4.2  Experimental  Data  and  Theoretical  Infotmation  Used  for  Parameter  Fitting. 

(1)  Structure.  The  structure  of  the  methoxy  radical  was  determined  by  microwave  spectral 
studies  of  Endo,  Saito,  and  Hirota;34  the  structural  parameters  are  given  in  Table  12.  Several  ab  initio 
calculations  have  been  done  that  predict  the  equilibrium  structure  of  methoxy  radical;37'40  however,  we 
have  used  the  experimental  values  for  the  equilibrium  structure.  The  structure  of  the  NO  radical  was 
taken  from  Herzberg41  and  is  also  given  in  Table  12. 
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Table  12.  Structural  Parameters  of  CH30  +  NO 


|  Internal  Coordinate1 

Expt’l.^ 

PES  1 

PES  2 

PES  3 

J  Methoxy  Radical  | 

R(CO) 

1.376 

1.376 

1.376 

1.376 

R(CH) 

1.0937 

1.080 

1.080 

1.080 

ZHCO 

108.25 

108.25 

108.25 

108.25 

ZHCH 

— 

110.66 

110.66 

110.66 

|  Nitrogen  Monoxide  Radical 

R(NO) 

1.1508 

1.1508 

1.1526 

1.1526 

*  Units  are:  R  in  A,  Z.  in  degrees. 
b  Methoxy  radical:  Endo,  Saito,  and  Hirota,  Ref.  36. 
e  Nitrogen  Monoxide  radical:  Herzberg,  Ref.  41. 


(2)  Vibrational  Analyses.  The  methoxy  radical  has  nine  vibrational  degrees  of  freedom. 
Assignments  of  the  fundamentals  from  infrared  spectral  studies  compiled  by  Jacox33  are  given  in 
Table  13.  The  fundamental  frequency  of  nitrogen  oxide  is  1,904.03  cm'1.41 

(3)  Bond  Dissociation  Energies.  As  with  the  product  fragments  CH3  and  N02>  the  bond 
dissociation  energies  have  not  been  firmly  established  for  the  methoxy  radical,  but  the 
thermochemistry  of  the  dissociation  of  methyl  nitrite  has  been  determined.21  Therefore,  we  have 
assumed  as  reasonable  bond  dissociation  energies  175.6  kcal/mol  and  61.0  kcal/mol  for  C-0  and  C-H 
bonds,  respectively,  and  assumed  a  bond  dissociation  energy  of  151  kcal/mol  for  N-0  to  give  the 
correct  thermochemistry. 

2.4.3  Results  of  the  Fit  The  internal  coordinates  of  the  methoxy  and  nitrogen  oxide  fragments 
predicted  by  these  potential  energy  functions  are  shown  for  comparison  with  the  experimental  values 
in  Table  12.  Assuming  the  bond  dissociation  energies  given  above,  all  other  parameters  were  fit  to 
give  agreement  with  the  spectroscopic  data.  Calculated  fundamental  frequencies  for  the  three 
potentials  are  given  in  Table  13  along  with  the  experimental  values.  The  equilibrium  geometries  and 
fundamental  vibrational  frequencies  do  not  differ  among  the  three  potential  energy  surfaces.  The 
calculated  frequencies  agree  well  with  the  experimental  values  except  for  the  5(CH3)  mode  at 
1,362  cm'1  (the  experimental  value).  The  three  potential  energy  surfaces  predict  a  higher  value  for  this 
mode  (1,431.8  cm'1).  The  best  fit  parameters  for  the  V,  are  given  in  Table  14. 
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Table  13.  Normal  Mode  Frequencies^  for  CH30  +  NO 


|  Mode 

ExpLb>6 

PES  1 

PES  2 

PES  3  S 

|  Methoxy  Radical  | 

1  «(CH) 

2840 

2873.6 

2873.6 

2873.6 

u(CH) 

2774 

2758.7 

2758.7 

2758.7 

5(HCH) 

1487 

1486.8 

1486.8 

1486.8 

5(CHj) 

1362 

1431.9 

1431.8 

1431.8 

v(CO) 

1047 

1005.3 

1005.3 

1005.3 

8(HCO) 

653 

654.7 

654.7 

654.7 

|  Nitrogen  Monoxide  Radical  | 

|  ^>(NO) 

1904.03 

1905.9 

1905.9 

1905.9 

*  Frequencies  are  in  cm'1. 
b  Methoxy  radical:  Jacox,  Ref.  33. 
e  Nitrogen  Monoxide  radical:  Herzberg,  Ref.  41. 


Table  14.  Potential  Energy  Parameters  for  CH,0  +  NO  Term 


Morse  Parameters 

Bond 

De(eV) 

a  (A1) 

R«,(A)  ! 

CH 

2.645232 

2.275 

CO 

7.614798 

1.705 

1.376 

NO 

6.548033 

2.76 

1.1508 

j  Harmonic  Oscillator  Parameters  || 

Angle 

k  (eV/rad2) 

^(degrees) 

ZHCH 

1.725 

110.66 

ZHCO 

0.805 

108.25 
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3.  FEATURES  OF  THE  POTENTIAL  ENERGY  SURFACE 


Table  15  gives  the  relative  energies  of  all  reactant  and  product  species  of  interest  in  this  study. 

In  addition  to  giving  reasonable  values  for  the  equilibrium  properties  of  all  stable  species  of  interest, 
the  three  surfaces  correctly  describe  the  thermochemistry  of  the  reactions.  The  values  of  the  relative 
energies  calculated  from  the  three  potential  energy  surfaces  described  above  agree  well  with 
experiment  and  with  each  other  except  for  the  transition  state.  Also,  the  often  dramatic  changes  in 
structure  of  the  species  upon  reaction  are  also  described.  For  example,  as  nitromethane  undergoes 
bond  scission  to  form  CH3  +  NO*,  the  methyl  group  "flattens"  from  a  tetrahedral  geometry  to  a  planar 
one  as  the  asymptote  is  approached.  Similarly,  the  ONO  angle  opens  up.  These  changes  are  apparent 
in  Figure  2,  which  illustrates  C-N  bond  scission  in  nitromethane.  For  each  point  generated  at  a  C-N 
bond  distance,  all  other  internal  coordinates  are  relaxed  to  an  equilibrium  position  (a  point  at  which 
the  gradient  in  these  coordinates  disappears).  Figures  2a,  2c,  and  2e  are  the  potential  energy  as  a 
function  of  the  C-N  bond  distance  in  nitromethane  calculated  from  PES  1,  PES  2,  and  PES  3, 
respectively.  Figures  2b,  2d,  and  2f  show  the  corresponding  change  in  HCH  and  ONO  angles  as  a 
function  of  R(C-N)  (solid  and  dashed  lines,  respectively).  In  agreement  with  experiment,2  there  is  no 
barrier  to  back  reaction  of  bond  scission  of  nitromethane.  The  parameters  in  SCHjN02  and  SCH3„N02 
were  adjusted  such  that  there  is  no  back  reaction  barrier  for  these  surfaces.  The  potential  energy 
curves  shown  in  Figures  2c  and  2e  are  not  as  smooth  as  the  curve  in  Figure  2a.  The  bumps  in  the 
curves  are  artifacts  due  to  the  choice  of  parameters  for  the  weighting  functions  for  PES  2  and  PES  3 
and  are  sufficiently  small  that  they  should  not  affect  the  dynamics  results. 

Figures  3  and  4  show  simple  C-0  bond  rupture  of  cis-  and  trans-methyl  nitrite,  respectively,  for 
the  three  potential  energy  surfaces.  As  in  Figure  2,  all  internal  coordinates  are  allowed  to  relax  to  an 
equilibrium  position  for  each  C-0  value  in  the  figures.  Figures  3a,  3c,  and  3e,  and  4a,  4c,  and  4e  are 
the  potential  energy  curves  vs.  R(C-O)  for  cis-  and  trans-methyl  nitrite,  respectively.  The  potential 
energy  curve  calculated  from  PES  1  (Figures  3a  and  4a)  have  no  barrier  to  back  reaction;  however,  the 
potential  curves  calculated  from  PES  2  and  PES  3  have  small  barriers.  The  barrier  heights  for  the 
back  reaction  for  cis  and  trans  conformers  on  PES  2  are  1.5  and  1.0  kcal/mol,  respectively  (Figures  3b 
and  4b).  The  barriers  to  back  reaction  for  the  cis-  and  trans-  conformers  on  PES  3  are  0.75  and  0.76 
kcal/mol,  respectively  (Figures  3c  and  4c).  This  is  a  result  of  the  choice  of  parameter  values  that  give 
the  desired  transition  state  energies.  Because  these  exit  channel  barriers  are  so  small,  and  enormous 
difficulty  was  encountered  in  obtaining  the  desired  transition  state  energies  for  PES  2  and  PES  3,  we 
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Table  15.  Relative  Energies*  of  Species  on  the  CHjN02  Potential  Energy  Surface 


Molecule 

PES  1 

PES  2 

PES  3 

Expt’l. 

CH3NO2 

0.0 

0.0 

0.0 

0.0 

cis-CHjONO 

2.30 

3.04 

3.04 

2.3b 

3.36 

3.36 

3.36 

3.1e 

CHjO  +  NO 

43.0 

43.0 

43.0 

43.5d 

CHj  +  NOj 

60.0 

60.0 

60.0 

59.5m 

T.S.* 

216.4 

55.1 

47.6 

55.5f 

*  Energies  in  kcal/mol. 

b  S.  W.  Benson,  Thermochemical  Kinetics  (Wiley,  New  York,  1976). 
c  Gwirm,  Anderson,  and  Stelman,  Ref.  26. 
d  Bait  and  Robinson,  Ref.  21. 

*  Transition  stale  for  isomerization  of  CH,N02  to  CH,ONO. 
f  Wodtke,  Hintsa,  and  Lee,  Ref.  2. 

opted  to  use  these  parameter  values.  Figures  3b,  3d,  3f,  4b,  4d,  and  4f  show  the  changes  in  the  HCH 
and  ONO  angles  (solid  and  dashed  lines,  respectively)  as  C-0  bond  rupture  occurs. 

Figures  5  and  6  show  simple  O-N  bond  rapture  of  cis-  and  trans-methyl  nitrite,  respectively. 
Although  all  three  surfaces  do  not  have  barriers  to  the  back  reaction,  PES  2  and  PES  3  show  a  sharp 
change  in  the  potential  energy  curve  at  approximately  1.7  A  which  is,  as  discussed  above,  attributed  to 
the  choice  for  parameter  values  used  to  lit  the  isomerization  barrier.  This  is  the  point  at  which  the 
contribution  of  the  methyl  nitrite  potential  energy  term  rapidly  approaches  zero,  and  the  contribution  of 
the  CHjO  +  NO  term  becomes  the  dominant  term  in  the  potential  energy  surface.  The  switch  from  the 
domination  of  the  methyl  nitrite  term  in  the  potential  to  that  of  the  CH30  +  NO  term  is  more  gradual 
and  takes  place  over  a  longer  O-N  bond  range  on  PES  1.  Figures  5b,  5d,  and  5f,  and  6b,  6d,  and  6f 
show  the  changes  in  the  HCH  and  HCO  angles  (solid  and  dashed  lines,  respectively)  as  the  O-N  bond 
becomes  large  for  cis-  and  trans-methyl  nitrite,  respectively.  Although  the  changes  in  these  geometries 
show  discontinuities,  note  that  the  changes  are  small  compared  to  those  for  the  formation  of 
CHj  +  NOj.  This  is  because  the  methyl  group  remains  in  a  tetrahedral  arrangement  as  the  methoxy 
radical  is  formed. 
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The  transition  state  for  isomerization  of  nitromethane  to  methyl  nitrite  has  been  located  on  each 
of  the  three  potential  energy  surfaces.  The  internal  coordinates  of  the  transition  state  on  each  potential 
energy  surface  are  given  in  Table  16.  The  energy  of  the  transition  state  on  PES  1  is  216.4  kcal/mol 
PES  2  has  a  transition  state  energy  of  S5.1  kcal/mol,  the  value  suggested  by  Wodtke  et  al.2  The 
energy  value  of  the  transition  state  on  PES  3  is  47.6  kcal/mol,  which  is  the  value  predicted  from  the 
MINDO/3  calculations  of  Dewar  and  Ritchie.16  The  geometries  of  the  transition  states  are  similar, 
except  that  the  CN  and  CO  bonds  are  from  0.3  to  0.7  A  longer  for  PES  2  and  PES  3  than  for  PES  1. 
The  frequencies  of  the  vibrational  modes  of  the  transition  states  are  given  in  Table  17;  a  substantially 
larger  frequency  for  the  transition  state  occurs  for  PES  1  compared  to  those  of  PES  2  and  PES  3. 

Figures  7,  8,  and  9  are  contour  plots  of  the  potential  energy  as  a  function  of  the  CNO  angle  and 
the  C-N  bond  distance  for  PES  1,  PES  2,  and  PES  3,  respectively.  For  each  grid  point  (i.e.,  C-N  bond 
distance  and  CNO  angle)  in  these  plots,  all  other  internal  coordinates  are  allowed  to  relax  to  the 
equilibrium  position.  Regions  which  energetically  describe  the  species  CHjN02,  CH3ONO,  and 
CH3  +  N02  are  all  shown  on  these  plots:  The  portion  of  the  figures  which  describes  nitromethane  is 
the  well  in  the  lower  right-hand  comer,  this  well  corresponds  to  the  equilibrium  C-N  bond  distance  in 
nitromethane  of  1.479  A  and  Bo*,  =  117.1°.  The  lowest  relative  energy  value  of  the  system  is  located 
at  this  grid  point,  and  has  the  value  of  0.0  kcal/mol  for  all  three  surfaces.  The  region  of  the  figure 
that  describes  CH3ONO  is  located  at  the  lower  left-hand  comer,  corresponding  to  a  C-N  bond  distance 
of  approximately  2.3  A  and  Sq©  =  32°,  the  equilibrium  C-N  bond  distance  and  CNO  angle  in  methyl 
nitrite.  The  value  of  the  lowest  energy  point  in  the  methyl  nitrite  well  is  3.4  kcal/mol  for  all  three 
surfaces,  indicating  that  this  is  the  trans  conformer.  The  flat  region  which  runs  across  the  upper 
portion  of  the  contour  plots  at  a  constant  value  of  60.0  kcal/mol  corresponds  to  the  CH3  +  N02  region 
of  the  surfaces. 

On  PES  1  (Figure  7),  the  barrier  to  isomerization  from  nitromethane  to  methyl  nitrite  is  located  at 
a  CNO  angle  of  approximately  75°  and  a  C-N  bond  distance  of  1.567  A,  and  has  a  value  of 
216.4  kcal/mol.  At  the  energies  considered  in  this  study  (190  kcal/mol  is  the  highest  total  energy  of 
the  system  and  includes  the  zero-point  energy  of  nitromethane),  isomerization  to  methyl  nitrite  from 
nitromethane  cannot  occur  on  PES  1  via  the  transition  state,  but  it  can  occur  if  a  trajectory  undergoes 
C-N  bond  scission  and  leaves  the  nitromethane  well,  passes  through  the  flat  region  of  the  potential  (the 
CHj  +  NO*  region),  and  then  moves  inward  to  the  methyl  nitrite  well.  For  this  surface,  isomerization 
is  a  two-step  process:  bond  scission  to  CH3  +  N02  followed  by  recombination  to  methyl  nitrite.  For 
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Table  16.  Structural  Parameters  of  the  Transition  State 


Internal  Coordinate* 

PES  1 

PES  2 

PES  3 

R(CN) 

1.567 

2.276 

2.196 

R(CO) 

1.713 

2.022 

1.950 

R(ON) 

1.236 

1.256 

1.260 

R(N=0) 

1.185 

1.189 

1.187 

R(CH)b 

1.078 

1.078 

1.078 

Z(ONO) 

119.19 

122.22 

122.45 

Z(CNO) 

74.30 

62.17 

61.90 

Z(CNO')c 

129.78 

131.25 

137.69 

Z(HoCiy4 

97.5 

107.7 

107.8 

98.8 

109.1 

109.3 

Z(H„CO)4 

114.9 

112.12 

112.02 

Z(H,CO) 

122.6 

118.8 

118.3 

x(CONO0 

122.3 

118.5 

118.0 

x(H,CON) 

70.46 

86.39 

87.70 

xd^CON)4 

127.18 

123.27 

130.84 

179.10 

178.86 

178.77 

92.70 

70.66 

69.74 

267.29 

288.66 

289.48 

*  Units  are:  R  in  A,  Z  in  degrees.  H,  designates  the  hydrogen  atom  in  the  plane  of  die  molecule,  and  H.  designates  the 
hydrogen  atom  out  of  the  plane  of  die  molecule. 

k  Does  not  differ  for  in-plane  and  out-of-plane  hydrogen  atoms. 

*  O'  denotes  die  terminal  oxygen. 

4  The  transition  state  stroctures  do  not  have  symmetry;  therefore,  the  H.CH,  and  H.CO  angles  differ  depending  on  relative 
orientation  to  the  CONO  frame.  Thus,  two  values  are  reported  for  both  angles  involving  the  two  H.  atoms. 
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Table  17.  Vibrational  Frequencies1  at  the  Transition  State 


PES  1 

PES  2 

PES  3 

2477.25 

3013.62 

3008.79 

2451.46 

3005.18 

2997.59 

2389.64 

2888.21 

2880.76 

1209.32 

1464.60 

1458.84 

1142.83 

1386.65 

1376.52 

1134.97 

1369.63 

1367.83 

1049.27 

1306.90 

1292.28 

1017.98 

1060.68 

1016.24 

720.33 

757.34 

730.74 

663.59 

541,17 

500.57 

605.40 

438.92 

378.74 

399.70 

405.33 

337.28 

320.39 

154.28 

140.80 

228.26 

140.30 

121.70 

3833.431 

1215.021 

1 127.571 

*  Frequencies  are  in  cm'1. 


33 


example,  suppose  a  molecule  undergoes  C-N  bond  scission  to  form  CH,  +  N02.  As  the  radicals 
recede  from  one  another  with  rotational  motion,  it  is  possible  that  one  of  the  oxygen  atoms  on  the 
rotating  NO}  can  come  within  interaction  range  of  the  carbon  atom  on  the  methyl  radical  during  this 
rotational  recession.  If  the  translational  energies  and  rotational  energies  of  the  receding  particles  are 
not  too  large,  the  ONO  group  will  recombine  with  the  methyl  group  via  one  of  the  oxygens  to  form 
methyl  nitrite.  Results  of  trajectories  run  on  PES  1  indicate  that  this  mechanism  is  reasonable;  they 
are  presented  in  the  Section  5. 

The  contour  plots  of  PES  2  and  PES  3,  shown  in  Figures  8  and  9,  respectively,  have  common 
features  except  for  the  barrier  heights.  The  transition  state  on  PES  2  (Figure  8)  has  a  value  of 
SS.l  kcalAnol,  and  is  located  at  a  C-N  distance  of  2.276  A  and  a  CNO  angle  of  62.2°.  The  transition 
state  on  PES  3,  shown  in  Figure  9,  is  located  at  a  C-N  bond  distance  of  2.196  A  and  a  CNO  angle  of 
61.9°,  and  has  an  energy  value  of  47.6  kcal/mol.  The  dissociation-recombination  mechanism  can  also 
occur  on  these  surfaces,  but  the  direct  isomerization  (i.e.,  concerted  C-N  bond  scission  and  C-0  bond 
formation  via  the  transition  state)  is  more  probable. 

4.  DETAILS  OF  THE  CALCULATIONS 

Ensembles  of  trajectories  with  initial  energies  in  the  range  64.6  to  161.4  kcalAnol  above  the  zero- 

a 

point  energy  of  nitromethane  (the  total  zero  point  energy  of  nitromethane  is  29.3  kcal/mol),  were 
integrated  using  a  variable-step  size  Adams-Moulton  fourth-order  predictor-corrector  integrator.42 
Relative  error  tolerance  was  set  at  10-7.  For  all  three  surfaces  and  for  most  energies,  Reactions  I,  II, 
and  III  were  observed. 

The  initial  conditions  were  selected  using  Metropolis  Monte  Carlo  sampling.43,44  Before  each 
ensemble  of  trajectories  is  integrated,  the  structure  of  the  nitromethane  molecule  is  set  at  the 
equilibrium  geometry,  and  the  desired  total  energy  of  the  system  in  the  form  of  kinetic  energy  is 
placed  in  the  normal  modes.  A  warm-up  Maikov  walk  of  1,000,000  steps  is  taken  to  randomize  the 
energy  of  the  system.  A  trajectory  is  then  integrated  until  the  end  tests  described  below  are  satisfied, 
or  until  the  trajectory  integration  exceeds  a  maximum  time  of  30  psec.  A  sequence  of  2,000  Maikov 
moves  is  taken  from  the  starting  point  of  the  previous  trajectory,  and  the  integration/Markov  walk 
pattern  is  repeated  until  1,000  trajectories  have  been  integrated.  Because  of  the  Monte  Carlo  sampling 
technique,  approximately  50%  of  the  1,000  trajectories  actually  integrated  were  recounted  and  included 
in  the  averaging.43'44 
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The  aid  tests  are  based  on  the  values  of  the  weighting  functions  described  in  Section  II.  If 
Schj+noj  >  Sms  and  all  other  S,  <  S^,  then  reaction  to  CH,  +  N02  has  occured  and  die  trajectory  is 
stopped.  If  ScHjQjsio  >  Sm„,  then  isomerization  is  noted,  but  the  trajectory  is  continued  until  additional 
end  tests  are  satisfied  or  the  maximum  integration  time  is  exceeded.  If  SCH30+N0  >  and  all  other  St 
<  Sife,  then  formation  of  CH30  +  NO  has  occurred  and  the  trajectory  is  stopped.  Because  the  value  of 
the  S)  are  dependent  on  geometric  factors,  these  end  tests  are  equivalent  to  using  internal  coordinates 
to  determine  when  reaction  has  occurred.  The  value  of  was  chosen  to  be  0.S  because  when  S,  = 
0.5,  then  the  contribution  of  the  corresponding  V,  to  the  total  potential  energy  is  approximately  50%. 

If  S,  >  0.5,  then  the  character  of  the  total  potential  energy  is  dominated  by  the  V,.  If  the  S,  =  0.0,  then 
the  V,  contributes  nothing  to  the  total  potential  energy.  Upon  decomposition  to  product  j,  all  St  (i*j) 
terms  approach  zero  and  Sj  approaches  unity  as  the  asymptotic  limit  of  the  receding  decomposition 
products  is  approached.  Rather  than  integrate  trajectories  until  this  limit  is  reached  and  all  V|S,  terms 
(i*j)  equal  zero,  we  found  that  when  all  S,  <  S^al  x  10"20,  reaction  had  definitely  occurred.  The  end 
tests  were  tested  extensively  by  running  the  same  ensemble  of  trajectories  for  a  full  30  psec  without 
invoking  end  tests  and  then  manually  determining  the  outcome  of  the  trajectories.  The  results  were 
the  same.  Additionally  tire  Markov  walk  is  restricted  such  that  Sq,^  must  have  a  value  greater  than 
0.01,  and  ScH30N0-  and  cannot  be  greater  than  S^.  These  restrictions  ensure  that  the  total 

potential  energy  at  the  beginning  of  each  trajectory  has  some  nitromethane  character  and  will  not  allow 
the  system  to  walk  into  the  methyl  nitrite  well  during  initial  condition  selection.  Sampling  of  phase 
space  in  this  manner,  however,  included  some  points  in  the  product  region  (CH,  +  NOj).  These  points 
corresponded  to  trajectories  which  did  not  undergo  at  least  one  C-N  vibration  as  nitromethane  before 
forming  CH,  +  NO,  and  were  not  included  in  the  ensemble  averages. 

5.  RESULTS  AND  DISCUSSION 

Before  presenting  and  discussing  the  results  of  the  calculations  of  ensembles  of  trajectories,  it  is 
worthwhile  to  consider  the  mechanisms  by  which  the  reactions  occur  on  the  three  potential  energy 
surfaces.  We  do  this  by  presenting  various  plots  of  selected  trajectories  that  are  typical  of  those  that 
give  the  ensemble  averaged  results  discussed  in  Section  5.2.  Two  types  of  plots  are  used  to  illustrate 
the  behavior  of  the  dissociation  and  isomerization  reactions.  The  first  type  are  plots  as  a  function  of 
time  of  the  weighting  functions  S|  that  are  involved  in  the  reaction.  These  weighting  functions,  which 
were  used  as  "switching"  functions  in  the  potential -energy  surface  formulations,  are  also  useful  for 
following  the  progress  of  reactions  in  the  trajectories.  The  other  type  are  plots  as  a  function  of  time  of 
the  internal  coordinates  that  are  involved  in  the  reaction. 
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5.1  Reaction  Mechanisms. 


5.12  PE S  1. 


CH«NO,  ->  CH,  ±  NO,  -»  CH,ONO  -»  CH,0  +  NO 

Figures  10a-d  and  lla-g  show  the  weighting  functions  and  internal  coordinates  for  a  typical 
trajectory  which  results  in  isomerization  to  methyl  nitrite  and  subsequent  dissociation  to  CH30  +  NO. 
The  total  energy  of  the  trajectory  is  69.2  kcal/mol  above  the  zero-point  energy  of  nitromethane 
(29.3  kcal/mol). 

Figures  lOa-d  are  the  weighting  functions  SCHjN02,  Schj+no^  Schjo-no.  31x1  schjo-*no.  respectively, 
plotted  as  functions  of  time.  For  the  first  0.6  psec  of  the  trajectory,  Sch^  is  unity,  while  all  other  S, 
are  zero  (with  a  few  minor  fluctuations).  sharply  decreases  to  zero  as  concurrently 

rises  to  one.  indicating  that  the  total  potential  energy  is  no  longer  characterized  mainly  by  the 
nitromethane  term  but  rather  by  the  terms  which  describe  methyl  and  nitrogen  dioxide  radicals.  After 
approximately  0.1  psec,  however,  drops  to  zero  as  rises  to  unity,  indicating  that  the 

geometry  of  the  system  is  consistent  with  that  of  methyl  nitrite  wherein  O'  is  bound  to  the  carbon 
atom  and  O  is  the  terminal  atom  in  methyl  nitrite.  At  this  time,  the  total  potential  energy  is  described 
by  Vchjo-no.  Schjo-no  is  unity  (with  minor  fluctuations)  for  approximately  1.75  psec,  after  which  it 
drops  to  zero  while  So^f  m  rises  to  one,  indicating  that  formation  of  CH30  +  NO  has  occurred. 
Therefore,  by  following  these  weighting  functions,  it  is  seen  that  the  reaction  mechanism  is: 

CH3N02  ->  CHj  +  NOj  -*  CHjONO  ->  CH30  +  NO. 

This  interpretation  of  the  reaction  is  confirmed  by  Figures  1  la-g  which  show  various  internal 
coordinates  vs.  time  for  the  same  trajectory.  Figure  1  la  shows  the  C-N  bond  distance  during  the 
trajectory.  For  the  first  0.6  psec,  the  C-N  bond  oscillates  about  its  equilibrium  nitromethane  value  of 
1.479  A.  It  then  begins  to  increase  (indicating  bond  scission)  but  then  decreases  and  oscillates  about 
2.3  A,  which  is  approximately  the  equilibrium  C-N  bond  distance  in  methyl  nitrite  until  the  end  of  the 
trajectory  (about  2.5  psec)  where  it  again  increases  until  it  meets  the  criterion  for  dissociation.  Figures 
lib  and  He  are  the  two  C-0  bond  distances  throughout  the  trajectory;  for  the  first  0.6  psec,  both  C-0 
bonds  oscillate  about  2.31  A,  the  equilibrium  C-0  bond  distance  in  nitromethane.  As  with  the  C-N 
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bond,  both  C-0  bond  distances  become  large  then  decrease  until  the  C-O'  bond  distance  oscillates 
about  1.421  A,  the  equilibrium  C-0  bond  distance  in  methyl  nitrite  (Figure  11c).  Figure  lid  shows 
the  ONO  angle  during  the  trajectory.  For  the  first  0.6  psec  of  the  reaction,  the  ONO  angle  oscillates 
about  the  nitromethane  equilibrium  angle  of  125.8°,  and  when  the  system  reaches  the  methyl  nitrite 
configuration,  it  oscillates  about  the  equilibrium  methyl  nitrite  value  (116.4°  or  118.65°,  depending  on 
whether  the  species  is  in  the  cis  or  trans  conformation).  Perhaps  the  most  dramatic  indication  of 
isomerization  is  the  CNO  angles  as  a  function  of  time,  shown  in  Figures  lle-f.  At  the  beginning  of 
the  trajectory,  both  CNO  and  CNO'  angles  oscillate  about  117.1°,  the  equilibrium  CNO  angle  in 
nitromethane.  After  isomerization  to  methyl  nitrite,  the  CNO'  angle  drops  from  117.1°  and  oscillates 
about  32°,  the  equilibrium  CNO  angle  in  methyl  nitrite  (Figure  Ilf).  Finally,  Figure  llg  shows  the 
00740  dihedral  angle  as  a  function  of  time  for  the  portion  of  the  trajectory  during  which  methyl 
nitrite  exists.  The  newly  formed  methyl  nitrite  undergoes  two  complete  rotations  and  then  settles  into 
the  cis  conformation  (0°)  until  dissociation  occurs. 

CHoNOo  ->  CH,  +  NO-  ->  CH,ONO  ->  CH,  +  NO, 

Figures  12a-c  and  13a-g  illustrate  a  trajectory  in  which  nitromethane  isomerizes  to  methyl  nitrite 
then  dissociates  to  CH3  +  N02.  The  total  energy  of  the  trajectory  is  115.3  kcal/mol  above  the  zero- 
point  energy  of  nitromethane.  Figures  12a-c  show  the  weighting  functions  SCH3N02,  SCH3+N02,  and 
Schjono-.  respectively;  all  other  S,  throughout  the  trajectory  are  zero.  For  the  first  2.5  psec  of  the 
trajectory,  the  SCHjNQ2  undergoes  large  oscillations  from  one  to  approximately  0.2;  concurrently, 

Sch3  +  no2  undergoes  large  oscillations  from  zero  to  0.8,  indicating  that  the  C-N  stretch  of  nitromethane 
is  highly  excited.  SCHjN02  then  drops  rapidly  to  zero  as  SCHj  +  ^  rises  to  unity,  indicating  C-N  bond 
scission.  After  another  0.1  psec,  +f(02  drops  to  zero,  and  SCHjONO-  rises  to  unity,  indicating  that 
isomerization  to  methyl  nitrite  has  occurred.  Large  oscillations  in  SCHj +NOj  and  SCH3ono-  31 
t  =  2.75  psec  indicate  that  the  C-0  stretch  is  highly  excited.  Dissociation  to  the  products  CH3+N02 
occurs  at  3.0  psec,  as  Schjono-  sharply  drops  to  zero  and  SCHj  +  NC>2  concurrently  rises  to  unity. 
Therefore,  the  reaction  can  be  summarized  as: 

CHjNOj  ->  CHj  +  N02  -»  CHjONO  -4  CH3  +  NO*. 

The  internal  coordinates  monitored  during  the  trajectory  are  illustrated  in  Figures  13a-g. 

Figure  13a  shows  the  C-N  bond  distance  throughout  the  trajectory  and  for  the  first  2.5  psec  of  the 
trajectory,  the  C-N  bond  is  highly  excited  as  indicated  by  the  large  oscillations  of  the  C-N  bond.  As 
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in  the  previous  illustration  (Figure  1  la),  the  C-N  bond  distance  becomes  very  large  at  2.5  psec,  then 
decreases  as  the  system  attains  the  methyl  nitrite  conformation.  Figures  13b-c  show  the  two  C-0 
bond  distances  as  a  function  of  time,  and  for  the  first  2.3  psec  of  the  trajectory  both  C-0  bonds 
oscillate  about  the  equilibrium  C-0  distance  in  nitromethane.  The  C-0  distances  undergo  large 
oscillations  and  then  become  very  large  at  the  same  time  the  C-N  distance  becomes  large,  but  then  the 
C-0  distances  decrease.  The  C-0  bond  oscillates  about  the  equilibrium  C-0  bond  distance  of  methyl 
nitrite  indicating  that  isomerization  has  occurred  while  the  C-O'  bond  oscillates  at  a  larger  value, 
indicating  that  O'  is  the  terminal  oxygen  of  the  methyl  nitrite  species.  Figure  13d  shows  the  ONO 
angle  during  the  trajectory.  At  the  beginning  of  the  trajectory,  the  ONO  angle  oscillates  about  125.8°, 
the  equilibrium  value  in  nitromethane,  but  upon  isomerization,  it  oscillates  about  the  equilibrium  value 
of  methyl  nitrite  (116.4°  or  118.65°,  depending  on  whether  the  species  is  in  the  cis  or  trans 
conformation).  As  in  the  previous  trajectory,  the  isomerization  is  most  dramatically  illustrated  by  the 
behavior  of  the  CNO  angles  during  the  trajectory.  Figures  13e-f  show  the  two  CNO  angles  in  time, 
and  upon  isomerization  the  CNO  angle  changes  its  oscillation  about  117.1°  (the  equilibrium  value  in 
nitromethane)  to  oscillations  about  32°  (the  equilibrium  value  in  methyl  nitrite).  Figure  13g  shows  the 
CONO'  dihedral  angle  as  a  function  of  time  for  the  portion  of  the  trajectory  for  which  methyl  nitrite 
exists.  The  newly  formed  methyl  nitrite  undergoes  one  full  rotation  before  dissociation  occurs. 

CFLNO-  — »  CH,  +  NO, 

For  completeness,  we  have  included  a  representative  trajectory  showing  the  simple  C-N  bond 
scission  in  nitromethane.  The  total  energy  of  the  system  is  138.4  kcal/mol  above  the  zero  point 
energy  of  nitromethane.  Figures  14a-b  are  plots  of  the  weighting  functions  SCH3N02  and  SCH3  +N02  as 
functions  of  time,  respectively;  all  other  S,  throughout  the  trajectory  are  zero  indicating  that  only  the 
terms  VCH3NQ2  and  +NOz  contribute  to  the  total  potential  energy  as  the  reaction  proceeds.  With  a 
few  small  oscillations,  SCH3NC>2  remains  unity  throughout  the  trajectory  and  then  drops  sharply  to  zero, 
while  Schj +N02  ^ses  to  unity  indicating  that  C-N  bond  scission  has  occurred.  Thus,  the  reaction  is; 

CHjN02  ->  CH3  +  N02. 

Figures  14c,  14d,  and  14e  show  the  C-N,  C-O,  and  C-O'  bond  distances,  respectively,  throughout 
the  trajectory;  each  bond  oscillates  about  its  equilibrium  nitromethane  value  until  the  bond  scission 
occurs,  at  which  time  the  C-N,  C-O,  and  C-O'  bond  distances  simultaneously  become  large.  Figures 
14f,  14g,  and  14h  show  the  ONO,  CNO,  and  CNO'  angles,  respectively,  during  this  trajectory.  The 
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angles  oscillate  about  their  equilibrium  values  in  nitromethane  (125.8°  for  ONO  and  117.1°  for  CNO) 
until  dissociation. 

The  trajectory  calculations  using  this  potential  energy  surface  indicate  that  isomerization  to 
methyl  nitrite  is  a  two-step  process: 

(1)  CHjN02  undergoes  C-N  bond  rupture  to  form  CH3  +  N02; 

(2)  CH3  and  N02  recombine  to  form  CH3ONO,  which  then  dissociates  to  either  CH30  +  NO  or 
CHj  +  N02. 


5.1.3  PES  2  and  PES  3.  Because  isomerization  on  both  of  these  surfaces  occurs  in  the  same 
manner,  we  have  chosen  to  illustrate  the  mechanism  of  isomerization  and  dissociation  using  PES  3. 
The  mechanism  of  isomerization  on  these  surfaces  is  concerted;  as  the  C-N  bond  breaks,  the  C-0  bond 
is  formed.  As  in  the  previous  section,  we  will  present  both  internal  coordinates  and  weighting 
functions  during  the  course  of  the  representative  trajectories. 

CH,NO,  ->  CH,ONO  -»  CH,0  +  NO 

Figures  15a-d  and  16a-g  show,  respectively,  the  weighting  functions  and  internal  coordinates 
during  a  trajectory  that  results  in  the  formation  of  CH30  +  NO.  The  total  energy  of  the  trajectory  is 
64.6  kcal/mol  above  the  zero  point  energy  of  nitromethane. 

For  the  first  6  psec,  SCHjN02  is  unity  with  minor  fluctuations,  indicating  that  the  total  potential  is 
described  by  the  nitromethane  term.  At  approximately  6  psec,  SCH307li0  shows  a  sharp  peak; 
immediately  thereafter  Sq^ono-  rises  to  unity  as  SCH3N02  correspondingly  goes  to  zero,  suggesting  that 
at  the  time  immediately  prior  to  isomerization  both  CNO'  and  CNO  bends  are  highly  excited.  Upon 
isomerization,  the  methyl  nitrite  species  remains  for  another  7  psec,  at  which  time  O-N  bond  scission 
occurs.  Therefore,  the  weighting  functions  indicate  that  the  mechanism  of  isomerization  and  formation 
of  CH30  +  NO  is: 


CHjN02  -»  CHjONO  ->  CH30  +  NO. 

Figures  16a,  16b,  and  16c  are  the  C-N,  C-O,  and  C-O'  distances  throughout  the  trajectory.  For 
the  first  1.75  psec,  the  C-N  and  both  C-O  bonds  are  oscillating  well  above  their  equilibrium  values  in 
nitromethane.  The  NQ  group  settles  down  for  about  2.0  psec,  oscillating  about  the  equilibrium  values 
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of  nitromethane.  The  C-N  and  C-0  bonds  then  undergo  large  oscillations  before  isomerization  at 
6.0  psec.  At  6.0  psec,  the  C-0  bond  begins  to  oscillate  about  the  equilibrium  bond  distance  in  methyl 
nitrite  (1.421  A)  (Figure  16b),  while  the  C-N  bond  oscillates  about  the  equilibrium  distances  in  methyl 
nitrite  and  the  C-O'  oscillates  at  a  large  bond  distance  indicating  the  O'  is  the  terminal  oxygen  of  the 
methyl  nitrite.  Figure  16d  shows  the  ONO  angle  during  the  trajectory;  for  the  first  6  psec  of  the 
trajectory,  the  ONO  angle  oscillates  about  its  equilibrium  value  in  methyl  nitrite.  Figures  16e  and  16f 
are  the  CNO  and  CNO'  angles  during  the  trajectory.  Figure  16e  is  the  most  dramatic  illustration  of 
isomerization:  The  CNO  angle  oscillates  about  117.1°  while  the  molecule  is  in  the  nitromethane 
configuration,  but  after  isomerization  it  oscillates  around  32°.  Figure  16g  is  the  CONO'  torsional 
angle  in  methyl  nitrite.  Upon  isomerization,  the  CONO'  undergoes  one  full  rotation,  settles  into  the 
cis  confoimer  for  approximately  1  psec,  undergoes  free  rotation  for  the  next  2  psec,  and  again  settles 
down  into  the  cis  conformation  before  once  again  undergoing  free  rotation.  At  12  psec,  die  methyl 
nitrite  converts  to  the  trans  conformation  before  finally  dissociating. 

CH,NO.  ->  CH,ONO  ->  CH,  +  NO, 

Figures  17a-d  and  18a-g  illustrate  the  weighting  functions  and  geometries  of  a  trajectory  which 
results  in  direct  isomerization  from  nitromethane  to  methyl  nitrite  and  subsequent  C-O  bond  scission  to 
form  CH3  +  N02.  Figures  17a-d  show  the  weighting  functions  SCHjN02,  SCH30N0.,  SCH301<o,  and  SCH,  +  N02, 
respectively.  The  total  energy  of  the  system  is  69.2  kcal/mol  above  the  zero  point  energy  of 
nitromethane.  For  the  first  4  psec  of  the  trajectory  (see  Figure  17a),  the  system  is  described  mainly  by 
the  nitromethane  terms  with  fluctuations  among  the  CH3ONO',  CHjO'NO  and  CH3  +  N02  potential 
terms  (Figures  17b,  17c,  and  17d,  respectively).  At  4  psec,  SCH3N02  decreases  to  zero  as  SCH30140  rises 
to  unity.  The  system  is  then  described  by  the  methyl  nitrite  term  for  the  next  5  psec,  wherein  a 
fluctuation  with  the  CH3  +  N02  term  occurs,  indicating  excitation  of  the  C-0  bond.  The  C-0  bond 
rupture  occurs  1  psec  later.  Therefore,  the  weighting  functions  show  the  mechanism  of  isomerization 
to  be 

CHjNOj  ->  CH3ONO  -»  CH3  +  N02. 

Figure  18  shows  the  internal  coordinates  during  the  trajectory.  Each  coordinate  shows  distinct 
changes  upon  isomerization.  As  in  tire  trajectories  discussed  above,  the  CN,  CO,  and  CO'  distances 
shown  in  Figures  18a,  18b,  and  18c,  respectively,  oscillate  about  their  equilibrium  nitromethane  values. 
At  4.0  psec,  the  C-O'  bond  in  methyl  nitrite  is  formed,  as  evidenced  by  the  C-O'  oscillation  about  the 
equilibrium  bond  distance  for  methyl  nitrite  (Figure  18c).  Figures  18d-f  are  the  ONO,  CNO,  and 
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CNO'  angles,  respectively.  All  angles  oscillate  about  their  nitromethane  equilibrium  values,  but  after 
isomerization  they  oscillate  about  the  methyl  nitrite  equilibrium  values.  The  CNO  angles  (see  Figure 
18e  and  18f)  show  dramatic  change  upon  isomerization.  The  CNO  angle  (see  Figure  18e)  shows  an 
interesting  pattern  of  oscillation  from  7-9  psec  that  is  due  to  the  torsional  motion.  Note  that  the 
torsion  angle  (Figure  18g)  during  that  time  period  is  highly  excited  and  "free"  internal  rotation  is 
occurring,  whereas  for  the  3  psec  before  and  1  psec  after  the  free  rotation,  the  methyl  nitrite  is  in  the 
trans  conformation. 

3.2.  Ensemble  Results.  Table  18  is  a  compendium  of  the  results  of  the  ensembles  of  trajectories 
integrated  at  different  energies  on  the  three  potential  energy  surfaces.  The  values  under  the  column 
labeled  "Rejected"  are  the  numbers  of  trajectories  that  were  not  included  in  the  averaging  because  they 
did  not  undergo  at  least  one  C-N  vibration  as  nitromethane  before  forming  CH3  +  N02.  The  number 
of  rejected  trajectories  increases  with  increasing  energy  on  all  three  potential  energy  surfaces.  We 
chose  to  sample  phase  space  that  included  some  points  in  the  product  region  (CH3  +  NOj)  and  then 
"reject"  those  points.  This  was  done  to  ensure  that  all  of  the  reactant  phase  space  was  covered  in  the 
sampling. 

The  curves  in  Figure  19  are  probabilities  of  the  occurrence  of  Reaction  III,  isomerization  to 
methyl  nitrite  followed  by  dissociation  to  CH3  +  NOj,  as  a  function  of  energy  for  the  three  potential- 
energy  surfaces.  For  the  most  part,  the  probability  of  Reaction  ni  increases  with  increasing  energy  of 
the  system  for  all  three  surfaces.  The  probability  of  Reaction  in  for  PES  3  (stars)  increases 
dramatically  from  64  to  100  kcal/mol  and  then  levels  off.  The  curve  which  describes  the  probability 
of  Reaction  in  for  PES  2  (crosses)  increases  monotonically  with  energy  up  to  150  kcal/mol  and  then 
drops  to  a  lower  value  at  the  highest  total  energy  studied  (161.4  kcal/mol  above  the  zero-point  energy 
of  nitromethane).  The  curve  for  PES  1  (circles)  shows  a  small  increase  in  the  probability  of 
Reaction  in  at  the  high  energy  end  of  the  curve.  Only  a  small  number  of  trajectories  result  in 
Reaction  in  on  PES  1,  and  thus  statistical  error  is  high.  This  is  expected  because  of  the  features  of 
PES  1  discussed  in  Section  ni. 

Figure  20  shows  the  probabilities  of  Reaction  II,  isomerization  to  methyl  nitrite  followed  by 
dissociation  to  CHjO  +  NO,  as  a  function  of  energy  for  the  three  potential  energy  surfaces.  The 
behavior  is  quite  different  among  the  three  surfaces.  The  probability  of  Reaction  II  on  PES  3  (stars) 
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Table  18.  Numbers  of  Trajectories  and  Reactions 
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decreases  with  increasing  energy  over  the  range  of  64  to  100  kcal/mol.  At  energies  above  100  kcal/mol, 
the  reaction  probability  appears  to  be  less  energy  dependent  The  probabilities  of  Reaction  II  on  PES 
2  and  PES  1  (denoted  by  crosses  and  circles,  respectively)  are  approximately  constant  over  the  energy 
range  studied. 

The  following  ratios  were  calculated: 

B,  =  Reaction  n/(Reaction  I  +  Reaction  III); 

Bj  =  (Reaction  n  +  Reaction  ni)/Reaction  I;  and 

B3  =  Reaction  II/Reaction  in. 

The  ratio  B,  corresponds  to  the  Wodtke  et  al.2  measurements  assuming  they  did  not  distinguish 
between  the  NO*  formed  in  Reaction  I  to  that  which  might  be  formed  in  Reaction  III.  However,  in 
estimating  the  isomerization  barrier  from  experimental  data,  it  seems  that  B2  would  be  the  ratio  which 
needs  to  be  considered.  B3  is  the  branching  ratio  of  the  decomposition  reactions  of  the  methyl  nitrite. 

Figure  21  shows  plots  of  the  branching  ratios  Bt  (circles)  and  B2  (crosses)  for  PES  1  (Figure  21a), 
PES  2,  (Figure  21b)  and  PES  3  (Figure  21c).  Figure  22  shows  plots  of  the  branching  ratio  B3  for  PES 
1  (Figure  22a),  PES  2  (Figure  22b)  and  PES  3  (Figure  22c).  The  Blt  Bj,  and  B3  for  PES  1  have  been 
included  for  completeness;  however,  the  statistical  errors  are  very  large  for  the  ratios  due  to  the  small 
number  of  isomerization  trajectories.  The  behavior  of  the  curves  reflect  the  energy  dependences  of  the 
reaction  probabilities  discussed  above.  However,  the  only  experimental  value  for  Bt  available  for 
comparison  corresponds  to  an  energy  near  the  threshold  of  C-N  bond  dissociation.  The  calculated 
value  at  this  energy  for  PES  2  is  0.032  whereas  the  value  for  PES  3  at  the  same  energy  is  0.79,  in 
reasonable  agreement  with  experiment  (0.60)2.  This  suggests  that  the  barrier  height  of  a  tight 
transition  state  such  as  those  described  by  PES  2  and  3  should  be  slightly  larger  than  47  kcal/mol. 

The  curves  for  B3  for  PES  2  and  PES  3  (Figure  22)  have  similar  energy  dependences.  At  lower 
energies  (below  100  kcal/mol),  B3  is  larger  than  1.0,  indicating  that  upon  isomerization,  O-N  bond 
scission  in  methyl  nitrite  is  more  probable  than  C-0  bond  scission.  However,  at  higher  energies,  B3  is 
approximately  constant  and  less  than  unity,  indicating  that  C-0  bond  scission  is  more  probable 
following  isomerization  than  O-N  bond  scissioa 
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Two  questions  arise  from  the  behavior  of  the  reaction  probabilities  and  the  branching  ratios: 

(1)  Why  does  the  C-0  bond  scission  (Reaction  ni)  occur  at  all  energies,  and  why  is  it  more 
probable  than  O-N  bond  scission  (Reaction  n)  in  methyl  nitrite  at  high  energies,  even  though  Reaction 
in  requires  16  kcal/mol  more  energy? 

(2)  Why  is  the  probability  of  Reaction  n  essentially  not  energy  dependent  on  PES  2  but 
decreases  as  a  function  of  energy  on  PES  3? 

Unimolecular  decomposition  of  methyl  nitrite  under  the  experimental  conditions  of  Wodtke  et  aL2 
would  most  likely  result  solely  in  O-N  bond  breaking  to  form  CH30  +  NO  because  of  the  large 
difference  in  the  dissociation  energies  of  the  O-N  and  C-0  bonds.  However,  the  C-0  bond  that  is 
formed  when  nitromethane  isomerizes  to  methyl  nitrite  is  very  likely  initially  highly  excited.  That  is, 
eneigy  is  selectively  deposited  in  this  bond  in  the  isomerization  reaction,  and  the  transfer  of  only  a 
small  amount  of  energy  from  other  molecular  modes  will  result  in  the  breaking  of  the  bond.  Thus, 
this  mode  selectivity  in  the  energy  distribution  of  the  nascent  methyl  nitrite  leads  to  the  rupture  of  the 
stronger  C-0  bond  (preferentially  at  high  energies)  seen  in  our  results. 

The  observed  increasing  probability  of  Reaction  III  with  increasing  energy  is  simply  due  to  die 
increasing  amount  of  energy  that  is  available  among  the  various  modes  for  transfer  to  the  C-0  bond. 
Why  then  is  Reaction  n  essentially  constant  as  a  function  of  energy  for  PES  2  but  not  for  PES  3? 

This  can  be  attributed  to  the  difference  in  the  barrier  heights  for  isomerization  from  nitromethane  to 
methyl  nitrite  for  the  two  surfaces.  The  barrier  height  for  isomerization  on  PES  2  is  53.1  kcal/mol  and 
47.6  kcal/mol  for  PES  3.  If  we  assume  that  all  of  the  energy  of  isomerization  goes  into  the  newly- 
formed  C-0  bond,  that  the  C-0  bond  on  PES  2  will  have  7.5  kcal/mol  more  energy  than  the  C-0 
bond  on  PES  3  because  of  the  difference  in  the  barrier  heights.  On  PES  2  only  2  kcal/mol  of  energy 
must  transfer  from  other  modes  to  provide  the  energy  needed  to  break  the  C-0  bond  (bond 
dissociation  energy  is  57  kcalAnol),  and  thus  the  reaction  is  essentially  independent  of  the  total  energy. 
However,  on  PES  3  9.4  kcal/mol  must  flow  into  the  C-0  bond  to  break  it  The  higher  energy 
requirement  for  C-0  bond  scission  on  PES  3  could  explain  the  observed  energy  dependence  of 
Reaction  II  on  this  surface. 

Figure  23  show  the  translational  energy  distribution  of  products  from  Reaction  I  (hollow  bars) 
and  Reaction  in  (filled  bars)  for  PES  3  at  69.2  (Figure  23a),  115.3  (Figure  23b),  and  161  kcal/mol 
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(Figure  23c)  above  the  zero-point  energy  of  nitromethane.  Although  there  are  no  exit  barriers  for 
Reaction  I  and  in  on  PES  3,  the  distributions  do  not  peak  at  zero  as  expected  for  simple  bond  scission 
reactions.  This  behavior  has  been  observed  in  classical  trajectory  studies  of  unimolecular  reactions  of 
NHj,  CH*.  and  SiH*.43 

There  is  an  increase  in  the  release  of  product  translational  energy  with  increasing  total  energy,  as 
shown  by  the  distributions  in  Figure  23  and  by  the  average  values  given  in  Table  19.  Note  that  the 
distribution  of  products  for  Reaction  III  are  strikingly  similar  to  those  of  Reaction  I. 

Figures  24  and  25  are  the  same  as  Figure  23,  except  they  are  for  the  product  translational  energy 
distributions  for  PES  2  and  PES  1,  respectively.  The  probability  of  the  formation  of  methyl  and 
nitrogen  dioxide  radicals  from  C-0  bond  dissociation  in  methyl  nitrite  (Reaction  III)  is  much  smaller 
than  that  of  Reaction  I  for  all  surfaces;  however,  the  product  translational  energy  distribution  for  both 
reactions  appear  to  have  the  same  features.  Figures  26  and  27  are  plots  of  the  average  relative 
translational  energies  of  the  products  for  Reactions  I,  II,  and  III  on  PES  3  and  PES  2,  respectively,  as 
functions  of  energy.  The  average  translational  energy  of  die  products  of  Reaction  II  are,  for  the  most 
part,  larger  than  the  average  translational  energy  of  the  products  of  Reaction  I  and  Reaction  III  at  all 
energies  studied.  This  is  expected  because  the  energy  required  for  Reaction  n  (O-N  bond  scission  in 
methyl  nitrite,  41  kcalAnol)  is  much  less  than  that  required  for  Reaction  I  or  Reaction  III 
(60  and  57  kcal/mol,  respectively).  The  molecules,  whether  nitromethane  or  methyl  nitrite,  have 
essentially  the  same  internal  energy  before  decomposition  and,  therefore,  products  of  Reaction  n  have 
at  least  17  kcal/mol  more  energy  to  be  distributed  among  internal  and  translational  degrees  of  freedom. 

Figures  28  and  29  are  plots  of  the  average  rotational  energy  of  the  methyl  radicals  formed  from 
Reaction  I  (circles)  and  Reaction  in  (crosses)  on  PES  3  and  PES  2,  respectively.  Figures  30  and  31 
are  plots  of  the  average  rotational  energy  of  the  nitrogen  dioxide  radicals  formed  from  Reaction  I  and 
Reaction  in  on  PES  3  and  PES  2,  respectively.  As  for  the  product  translational  energy  distributions, 
remarkable  similarities  are  apparent  in  the  average  translational  and  rotational  energies  of  the  product 
of  Reactions  I  and  m  shown  in  Figures  26-31.  These  similarities  indicate  that  products  from 
Reactions  I  and  HI  might  be  very  difficult  to  distinguish  experimentally. 

Figure  32  shows  the  product  translational  energy  distribution  of  products  for  Reaction  II  for  PES  2 
(Figure  32a)  and  PES  3  (figure  32b)  at  69.2  kcalAnol.  The  distributions  are  quite  dissimilar,  however, 
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Table  19.  Average  Product  Translational  Energies 
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this  could  be  due  to  the  small  number  of  trajectories  resulting  in  Reaction  II  on  PES  2  at  this  energy, 
thus  giving  large  statistical  error.  The  distribution  in  Figure  32b,  however,  has  a  shape  very  similar  to 
the  translational  energy  distributions  of  the  products  of  Reactions  I  and  m  (Figures  23-25).  Tables 
20  and  21  give  the  average  rotational  energies  of  products  from  Reactions  I  and  m  and  Reaction  II, 
respectively.  The  internal  energy  of  all  decomposition  products  increase  with  increasing  energy. 
Exceptions  appear  in  Reactions  II  and  m  on  PES  1  and  are  attributed  to  the  poor  statistics. 

6.  SUMMARY 

We  have  performed  a  classical  dynamics  study  of  the  unimolecular  decomposition  of 
nitromethane.  Three  model  potential  energy  surfaces  were  used  in  the  study.  The  surfaces  differ 
mainly  in  the  heights  of  the  barrier  to  the  isomerization  of  nitromethane  to  methyl  nitrite.  PES  1  has  a 
very  high  barrier  to  isomerization  (216.4  kcal/mol)  whereas  the  barrier  for  PES  2  is 
55.1  kcal/mol,  the  value  suggested  by  Wodtke,  Hints  a,  and  Lee2  and  PES  3  has  a  barrier  of 
47.6  kcal/mol,  which  compares  well  with  the  predicted  MINDO/3  value  of  Dewar  and  Ritchie.16  Two 
primary  unimolecular  decomposition  pathways  were  suggested  by  Wodtke  et  al.2  in  interpreting  the 
results  of  their  infrared  multiphoton  dissociation  experiments: 


[I]  CH3N02  ->  CHJfOl  -*  CH3  +  N02  ; 

[H]  CH3  NO\  ->  CHtONO'  ->  CH20  ♦  NO  . 

An  additional  primary  unimolecular  decomposition  pathway  was  observed  in  the  trajectory  results  for 
all  three  potential  energy  surfaces, 

[nn  chjioI  ->  ch3ono '  ch3  *  no2  . 


Reaction  I  occurs  by  the  same  mechanism  for  all  three  potential  energy  surfaces.  The  isomerizing 
reactions,  Reactions  II  and  ID,  however,  occur  by  different  mechanisms  for  PES  1  and  PES  2  and  3. 
The  energy  of  the  barrier  to  isomerization  on  PES  1  is  larger  than  any  of  the  total  energies  of  the 
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Table  21.  Average  Rotational  Energies*  of  Products  CH30  +  NO 
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Energies  are  in  kcaJAnol  and  do  not  include  the  total  zero  point  energy  of  ni  tram  ethane  (293  kcal/mol). 
Trajectories  were  not  calculated  at  this  energy  on  PES  1. 


trajectories  calculated,  precluding  isomerization  via  the  transition  state  region  on  this  surface. 

However,  isomerization  occurs  on  this  surface  by  a  mechanism  of  dissociation  of  nitrom ethane  to 
CH3  +  N02,  followed  by  recombination  of  the  radicals  to  form  methyl  nitrite.  The  probability  of 
isomerization  on  this  surface  increases  with  increasing  energy. 

The  mechanism  of  isomerization  for  PES  2  and  PES  3  is  a  concerted  C-N  bond-breaking  and 
C-0  bond-forming  process.  PES  2,  the  surface  which  has  the  barrier  height  suggested  by  Wodtke  et  al.2, 
predicts  a  branching  ratio  at  least  an  order  of  magnitude  smaller  than  that  determined  from  the 
experimental  results.  PES  3,  the  surface  with  the  lowest  barrier  to  isomerization,  predicts  a  branching 
ratio  in  reasonable  agreement  with  that  reported  by  Wodtke  et  aL2 

Although  PES  2  and  PES  3  differ  by  only  a  few  kcal/mol  in  barrier  heights  and  show  the  same 
mechanisms  for  isomerization  and  decomposition,  the  dynamics  results  are  dramatically  different  The 
energy-dependent  branching  ratios  of  isomerization  reactions  to  C-N  bond  scission  reactions  for  PES  2 
increase  with  increasing  energy.  The  same  branching  ratios  for  PES  3  decrease  with  increasing 
energies  from  64.1  kcal/mol  to  103  kcal/mol,  then  remain  constant  up  to  161.4  kcal/mol.  There  is  a 
strong  energy  dependence  of  the  probability  of  Reaction  II  for  PES  3,  whereas  the  probability  of 
Reaction  II  remains  approximately  constant  as  a  function  of  energy  on  PES  2.  However,  the 
probability  of  Reaction  III  increases  with  increasing  energy  for  both  PES  2  and  PES  3. 

Average  translational  energies  of  the  products  for  Reaction  n  are  higher  than  those  of  Reactions  I 
and  III  at  all  energies.  The  product  translational  energy  distributions,  as  well  as  the  average 
translational  and  rotational  energies  of  the  products  of  Reactions  I  and  III  are  similar  for  all  three 
surfaces,  indicating  that  the  products  of  Reaction  I  and  Reaction  III  might  be  indistinguishable  in 
experiments. 

Although  the  trajectory  results  suggest  that  the  barrier  height  of  the  transition  state  should  be 
slightly  higher  than  47  kcal/mol,  several  good  quality  ab  initio  calculations14,15,17  predict  a  much  higher 
barrier  to  isomerization  for  a  tight  transition  state  of  the  type  described  by  the  potential  energy 
surfaces  in  this  study.  It  is  not  known  whether  a  much  better  electronic  structure  calculation  would 
lower  the  energy  value  of  this  transition  state  by  20  to  30  kcal/mol.  However,  there  have  been 
suggestions  that  another  transition  state  might  exist  which  resembles  a  van  der  Waals  complex.1*19 
Isomerization  through  this  transition  state  would  resemble  the  dissociation-recombination  mechanism 
seen  in  the  results  for  PES  1. 
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Although  gains  have  been  made  previously  in  understanding  the  electronic  structure  and  dynamics 
of  this  system,  this  study  uncovers  complexities  in  the  unimolecular  decomposition  of  nitromethane. 
Our  models  predict  that  nitromethane  has  three  decomposition  pathways  rather  than  the  two  observed 
in  the  experiment.2  Therefore,  experiments  need  to  be  performed  on  this  system  in  an  effort  to 
determine  the  occurrence  of  Reaction  m.  If  Reaction  in  is  found  to  exist,  this  will  affect  the  analysis 
used  to  estimate  the  isomerization  barrier.  Additionally,  our  models  indicate  that  there  are  two 
possible  mechanisms  for  isomerization  of  nitromethane  to  methyl  nitrite.  An  accurate  description  of 
the  transition  state  from  higher  quality  ab  initio  quantum  chemical  calculations  could  help  determine 
which  isomerization  mechanism  is  most  probable. 
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Figure  2.  Tod  Half:  Potential  Energy  as  a  Function  of  the  C-N  Bond  Distance  in  Nitromcthane  for  (a)  PES  1.  (c)  PES  2,  and  (e)  PES  3. 
Other  Internal  Coordinates  Are  Optimized  to  Minimize  the  Energy  at  Each  C-N  Distance.  Bottom  Half:  Optimized  HCH  angl 
(Solid  Line)  and  ONO  Angle  (Dashed  Line)  as  a  Function  of  C-N  Bond  Distance  for  (bl  PES  1.  (d)  PES  2.  and  (f)  PES  3. 


100- 


Figuie  3.  Too  Half:  Potential  Energy  as  a  Function  of  the  C-0  Bond  Distance  in  Cis-Mcthvl  Nitrite  for  (a)  PES  1,  (c)  PES  2.  and  (e)  PES  3. 
All  Other  Internal  Coordinates  Are  Optimized  to  Minimize  the  Energy  at  Each  C-0  Distance.  Bottom  Half:  Optimized  HCH 
Angle  (Solid  Line)  and  ONO  Angle  (Dashed  Line)  as  a  Function  of  C-O  Bond  Distance  for  (b)  PES  1,  (J)  PES  2.  and  (0  PES  3. 
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fe)  PES  3.  All  Other  Internal  Coordinates  Are  Optimized  to  Minimize  the  Energy  at  Each  0-N  Distance.  Bottom  Half: 
Optimized  HCH  Angie  (Solid  Line)  and  HCO  Angle  (Dashed  Line)  as  a  Function  of  0-N  Bond  Distance  for  (b)  PES  1. 
d)  PES  2.  and  (Q  PES  3. 
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Igure  6.  Top  Half:  Potential  Energy  as  a  Function  of  the  0-N  Bond  Distance  in  Trans-Methyl  Nitrite  For  (a)  PES  1.  (c)  PES  2. 
(e)  PES  3.  All  Other  Internal  Coordinates  Are  Optimized  to  Minimize  the  Energy  at  Each  0-N  Distance.  Bottom  Half: 
Optimized  HCH  Angle  (Solid  Line)  and  HCO  Angle  (Dashed  Line)  as  a  Function  of  0-N  Bond  Distance  for  (b)  PES  1 
(d)  PES  2.  and  (0  PES  3. 


Contour  Plot  of  Potential  Energy  (PES  2)  as  a  Function  of  the  C-N  Bond  Distance  and  CNO  Angle.  All  Other  Internal  Coordinates 
Are  Optimized  to  Minimize  the  Energy  at  Each  C-N.  CNO  Grid  Point.  Energy  Contour  Values  Are  in  kcal/mol.  and  Contour 
Values  From  0.0  to  75  kcal/mol  at  5  kcal/mol  Increments  Are  Shown. 


Are  Optimized  lo  Minimize  the  Energy  at  Each  C-N.  CNO  Grid  Point.  Energy  Contour  Values  are  in  kcal/mol.  and  Contour 
Values  From  0.0  to  75  kcalAnol  at  5  keal/mol  Increments  Are  Shown. 
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Figure  13. 
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Figure  15.  WcjfihtinR  Function?  (a)  SrH^0,iibISrH,n^,  (c)  Sru^n-  and  (d)  During  a 

Trajectory  Integrated  on  PES  3  Which  Results  in  Formation  of  Methvl  Nitrite  and 
Subsequent  Dissociation  to  CH,Q  +  NO  Via  a  Direct  Isomerization  Mechanism. 
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Figure  16. 
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Distance:  (b)  C-0  Bond  Distance:  (c)  C-0'  Bond  Distance:  (d)  ONO  Angle:  (e)  CNO 
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ZCNO  2  LC NO  1  /ONO 


PROBABILITY 


Figure  19.  Probability  of  Reaction  III  Occurring  as  a  Function  of  Energy  for  PES  1  (Circles).  PES  2 
(Crosses),  and  PES  3  (Stars).  The  Energy  Alone  the  Abscissa  Is  the  Total  Energy  Above 
the  Zero  Point  Energy  of  Nitromethane. 
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Figure  20.  Probability  of  Reaction  II  occurring  as  a  Function  of  Enerev  for  PES  1  (Circles).  PES  2 
(Crosses),  and  PES  3  (Stars).  The  Energy  Alone  the  Abscissa  Is  the  Total  Energy  Abov 
the  Zero  Point  Energy  of  Nitromethane. 
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Figure  27.  Average  Relative  Translational  Energy  of  the  Products  From  Reactions  I  (Circles).  II 
(Crosses),  and  m  (Stars)  on  PES  2  as  a  Function  of  Energy.  The  Energy  Along  the 
Abscissa  Is  the  Total  Energy  Above  the  Zero  Point  Energy  of  Nitromethane. 
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Figure  28.  Average  Rotational  Energy  of  Methyl  Radical  Formed  From  Reactions  I  (Circles)  and  III 
(Crosses')  on  PES  3  as  a  Function  of  Energy.  The  Energy  Along  the  Abscissa  Is  the 
Total  Energy  Above  the  Zero  Point  Energy  of  Nitromethane. 
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Figure  29.  Average  Rotational  Energy  of  Methyl  Radical  Formed  From  Reactions  I  (Circles)  and  III 
(Crosses)  on  PES  2  as  a  Function  of  Energy.  The  Energy  Along  the  Abscissa  Is  the 
Total  Energy  Above  the  Zero  Point  Energy  of  Nitromethane. 
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<E(ROT)>  (kcal/mol) 


Figure  30.  Average  Rotational  Energy  of  Nitrogen  Dioxide  Radical  Formed  From  Reaction  I  (Circles) 
and  III  (Crosses)  on  PES  3  as  a  Function  of  Energy.  The  Energy  Along  the  Abscissa  Is 
the  Total  Energy  Above  the  Zero  Point  Energy  of  Nitromethane. 
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Figure  31.  Average  Rotational  Fnerev  of  Nitrogen  Dioxide  Radical  Formed  From  Reaction  I  (Circles' 
and  III  (Crosses)  on  PES  2  as  a  Function  of  Energy.  The  Energy  Along  the  Abscissa  Is 
the  Total  Energy  Above  the  Zero  Point  Enenrv  of  Nitromethane. 
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Figure  32.  Distribution  of  the  Product  Relative  Translational  Energies  for  CH,0  +  NO  Formed  From 
Reaction  II  on  (a)  PES  2  and  (b)  PES  3  at  69.2  kcal/mol  Above  the  Zero  Point  Energy  of 

Nitromethane. 
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